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PREFACE 


In  April  1975>  the  Staff  Meteorology  Office  of  the  Air  Force  Weapons  Laborato)  r 
(AFWL/WE) , working  in  conjunction  with  the  Laser  Systems  Analysis  Branch  (AFWL/POA,, 
began  a special  project,  "Weather  Studies."  The  objective  of  this  project  was  to 
determine  the  potential  influence  of  the  atmosphere  and  its  associated  phenomena 
various  airborne*  applications  of  laser  devices.  As  the  study  progressed,  the  nut  i >rs 
became  Increasingly  Interested  in  the  problem  of  predicting  the  probability  of  a 
cloud-free  line-of-sight  (CFLOS)  between  two  moving  points.  Although  techniques 
were  available  for  handling  CFLOS  questions  of  a static  or  instantaneous  nature,  U 
became  apparent  that  these  techniques  were  not  designed  to  answer  the  kinds  of 
questions  being  raised  in  "Weather  Studies"  and  that  new  techniques  were  requir'd 

The  computer  model  presented  herein  is  a preliminary  step  in  developing  the  new 
techniques  required  to  solve  CFLOS  questions  involving  a moving  observer,  movir 
observed  point,  and/or  a time  dimension.  The  major  part  of  this  investigation  was 
performed  by  Ronald  J.  Kelson.  The  major  part  of  the  computer  model  was  perfo  mod 
by  Mead  B.  Wetherbc. 

The  authors  gratefully  acknowledge  the  constructive  criticisms  and  support  of 
their  co-workers  at  AFWL.  Thanks  are  due  especially  to  Captains  Steven  Edclman, 

John  King,  and  Carl  Curatola  for  their  valuable  contributions  during  the  development 
of  the  studies  and  analyses  presented  In  this  report. 
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SOME  ASPECTS  OF  ESTIMATING  THE  PROBABILITY 
OF  CLOUD- FREE  LINES-OF-SIGHT  TN  DYNAMIC  SITUATIONS 


SECTION  A — INTRODUCTION 

With  the  Introduction  of  optical  and  ini" rare ..  detection  and  tracking  devices  into 
the  Department  of  Defense  inventory,  interest  grew  in  developing  methods  for  esti- 
mating the  usability  of  such  devices  in  the  presence  of  clouds.  One  result  of  this 
interest  was  the  development  of  various  methods  for  predicting  the  probability  of 
having  a cloud-free  line-of-sight  (CFLOS)  between  two  points.  These  methods,  which 
have  been  available  in  various  forms  since  the  1960s,  provide  tools  for  estimating 
CFLOS  probabilities  between  a stationary  observer  and  a stationary  observed  point  at 
a specified  elevation  angle  in  the  presence  of  a given  amount  of  intervening  clouds. 
None  of  these  methods  (e.gv,  Lund  and  Shanklln  [1],  McCabe  [2],  ar.d  Rapp,  et  al  [3]) 
was  designed  to  address  CFLOS  questions  when  moving  observers,  moving  observed 
points,  or  time  dimensions  wore  involved.  These  kinds  of  questions  are  "dynamic 
CFLOS  questions." 

Many  USAF  applications,  actual  and  potential,  involve  the  use  of  optical  or 
Infrared  systems  or.  aircraft.  Some  applications  might  involve  an  airborne  system, 
a.g.,  as  a detector  of  a stationary  target  on  the  ground.  In  such  a case  one  might 
wish  to  predict  the  probability  that  chore  will  be  a CFLOS  between  the  detector  and 
the  target.  A typical  approach  to  this  problem  is  to  specify  the  altitude  at  which, 
and  the  range  of  depression  angles  within  which,  the  system  might  operate.  Then, 
given  the  climatological  frequency  distribution  of  cloud  conditions  for  the  area 
over  which  the  system  is  to  ce  employed,  use  some  technique,  such  as  that  outlined 
in  Rapp,  et  al.,  (3]>  to  predict  the  probability  of  having  a CFLOS  between  the  two 
vehicles  (points).  Since  at  any  instant  the  two  points  may  be  considered  to  be 
stationary,  the  approach  seems  reasonable  for  predicting  the  probability  of  an  instan 
tnneous  CFLOS.  Suppose,  however,  that  one  is  interested  not  in  some  "instantaneous" 
CFLOS,  but  rather  in  a CFLOS  which  lasts  for  some  finite  time,  on  the  order  of  one 
second.  Or,  suppose  that  both  the  observer  and  target  are  moving  and  that  one  wishes 
to  predict  the  probability  of  a CFLOS  lasting  for  some  finite  time  between  the  two. 
The  approach  outlined  above  cannot  handle  these  dynamic  problems. 

The  remainder  of  this  report  presents  the  results  of  the  work  to  date  on  dynamic 
CFLOS  problems  and  makes  available  the  computer  model  which  we  have  used  to  gain 
insight  into  problems  of  this  variety.  It  is  by  no  means  an  exhaustive  dissertation 
but  the  authors  hope  it  can  serve  as  a departure  point  for  others  .interested  in 
dynamic  CFLOS  problems. 
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SECTION  B — DYNAMIC  CFLOS  PROBLEMS 


Discussion  of  the  Variables 

One  of  the  mart  Important  "variables"  involved  In  dynamic  CFLOS  problems  is  the 
statement  of  the  problem.  The  following  problem  is  used  as  a basis  for  discussing 
the  variables  associated  with  dynamic  CFIDS  problems  in  general: 

Assume  a vehicle  B will  be  flying  at  a specific  altitude  and 
velocity  over  a specified  area.  A second,  faster  vehicle  A will 
be  introduced  at  an  arbitrary  initial  slant  range,  azimuth,  and 
elevation  angle.  The  second  vehicle  will  move  on  a straight- 
line  course  (Direct  Intercept  Course)  designed  to  effect  a 
collision  with  the  first  vehicle  Given  sufficiently  detailed 
cloud  in format ton  for  the  specified  area,  predict  the  probability 
that,  at  any  time,  a CFLOS  will  exist  between  the  two  vehicles 
for  at  least  t units  of  time  before  the  two  vehicles  reach  a 
given  slant- range  separation. 

The  variables  of  the  problem  are  ns  follows: 


a.  The  initial  positions  of  the  two  vehicles. 

b.  The  relative  speeds  and  paths  of  the  vehicles. 

c.  The  "boundary  conditions"  for  the  time  duration  ( ct ) of  the  CFLOS  and  the 
final  separation  distance  before  which  the  specified  CFLOS  must  occur. 

d.  The  cloud  characteristics  for  the  area  of  interest. 

Examination  of  the  Relative  Speeds  and  Paths  * # 

The  first  step  in  solving  dynamic  CFLOS  problems  is  to  determine  the  paths  to  be 
followed  by  the  vehicles.  Figure  1 depicts  the  geometry  of  the  situation  described 
above. 


A . B_:  Initial  positions  of  A and  B 
0 0 

(*A5  °ir  initial  azimuths  of  A and  B 

D : Initial  slant  range  separating 
A and  B 

D. : Distance  A must  travel  to  meet 
B 

a : Elevation  angle  from  B to  A 


f 

Figure  1.  Geometry  of  the  Direct  Inter- 
cept Course. 
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IT  VA,  \'B  represent  the  speeds  of  A and  B,  respectively,  then  the  relationship 
of  Da,  the  distance  A must  travel  to  meet  d,  to  the  other  parameters  (see  Appendix  A 
for  details)  is: 


whe  re 


and 


DA  » D • t>V> 


A 


(1) 


I?) 


F^  » Cos  (Cos  ('A  Sin  *>A  Tan 


»V)  Cos  I'., 


(o) 


P's  total  travel  distance  to  the  collision  point  (D.,)  is  directly  proportional 
to  A's  and  to  the  speed  ratios  involved.  That  is, 


Va 


({>) 


Knowing  the  total  distance  Involved,  one  can  use  direction  cosines  to  determine 
the  paired  x,y,c  positions  of  A and  B in  whatever  time  increments  are  desired.  The 
paths,  therefore,  are  defined. 

The  Boundary  Conditions  for  t and  S 

It  would  seem  that  one  should  now  be  able  to  specify  the  boundary  conditions  for 
the  time  duration  of  the  OFLOS  and  the  final  separation  distance  of  interest  and 
using  the  detailed  cloud  information  for  the  specified  area,  proceed  with  the  CFLOB 
analysis.  In  the  early  stages  of  the  study  the  authors  proceeded  in  Just  that  way, 
utilising  a randomly  generated  cloud  field  in  the  computer  model  as  a substitute  for 
the  "detailed  cloud  Information." 

Fortunately,  the  program  was  written  so  that  one  could  obtain  a highly  detailed 
output  of  the  progress  of  the  analysis.  For  a given  set  of  cloud  and  boundary  con- 
ditions, If  B was  repeatedly  reset  to  its  starting  position  and  A was  Introduced  at 
various  asimuths,  the  ratio  of  free/obstructed  lines-of-s.lght  was  almost  always 
greater  for  approaches  which  came  toward  B from  behind.  The  reason  for  this  result 
will  become  apparent  through  an  examination  of  the  concept  of  a "critical  cloud 
base.” 
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Concept  of  a Critical  Cloud  Base 

Figure  2 shows  the  geometry  to  be  considered  for  this  concept.  A and  D will 
again  travel  along  straight  lines  with  A following  the  direct  intercept  course. 


A , B : Initial  positions 
of  A and  B 

D : Initial  slant 

range  separating 
A and  B 

a : Elevation  angle 
from  B to  A 

r’  : Relative  azimuth 

between  A„  and  B . 

o oJ 

- PB!> 

Z : Cloud  base  height 


Figure  2.  Geometry  of  the  Critical  Cloud  Base. 


Consider  the  following  problem: 

For  the  initial  conditions  specified  (i-e  , VA/VB>  l),  a,  etc.), 
determine  that  height  above  the  plane  ox''  B,  (Z^),  above  which  the 
existence  of  clouds  w.ll  hove  no  effect  on  the  probability  of  having 
a C PIi)S  for  At  units  of  time  at  least  once  before  A and  B approach 
within  S of  each  other. 

The  detailed  development  given  in  Appendix  B shows  that 


z zv 


Sin 


a (s 


S + At 


YA\ 

¥~j 

*2 


(3) 


where  Zc  is  a critical  relative  cloud  height  ,».n  that  if  all  clouds  are  based  at 
Z > + 2^,  the  required  CFIDS  will  occur  &t  least  once.  Equation  (B)  shows  that 

one  should  not  specify  cloud  characteristics  for  any  model  situation  without  con- 
sidering Uv  "boundary  conditions"  and  dynamics  of  the  situation. 
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Figure  It.  Variation  of  Zc  with  Relative  Azimuth. 

travelling  at  a speed  ratio  (VA/VB)  of  3 to  1,  for  all  elevation  angles  less  than 
about  24°  at.  any  relative  azimuth.  Similarly,  2000-meter  cloud  bases,  for  the  con- 
ditions specified  on  Figure  3(d)  will  Insure  the  required  CFLOS  for  an  elevation 
angle  of  70°  for  relative  azimuths  of  approximately  109°  * P * l80°.  Note  that  the 
0-180°  relative  azimuths  can  be  "mirrored"  into  the  l80°-260°  azimuth  range  and  that 
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"elevation  angle"  could  as  easily  be  interpreted  as  depression  angle.  As  a result, 
can  be  treated  as  a cloud  separation  distance  and  therefore  could  be  applied  to 
cloud  tops. 

Figures  4(a)  through  4(d)  depict  the  variation  of  Zc  as  a function  of  initial 
relative  azimuth  and  elevation  angle  for  various  boundary  conditions.  The  value  of 
Zc  generally  increases  with  increasing  elevation  angle  and  .s  higher  for  frontal 
than  for  rearward  approaches.  Note,  however,  that  this  generalization  breaks  down 
when  the  approach  combines  large  elevation  angles  and  small  relative  azimuths. 

These  variations  should  have  a direct  bearing  on  how  one  compiles  the  statistics 
in  any  modeled  situation  for  a given  cloud  base.  If,  for  example,  the  CFLOS  prob- 
ability statistics  associated  with  a given  elevation  angle  were  grouped  into  a 
single  value  for  all  azimuths,  that  value  would  implicitly  be  weighted  by  the  higher 
values  associated  with  the  rearward  approaches.  On  the  other  hand,  grouping  the 
statistics  at  a given  azimuth  for  all  elevation  angles  would  implicitly  weight  the 
answer  by  the  higher  values  associated  with  the  lower  elevation  angles.  The  physi- 
cal reasons  are  apparent: 

a.  For  a given  elevation  angle  the  closure  rate  is  lower  for  large  relative 
azimuths. 

b.  For  a give:,  azimuth  the  closure  rate  is  lower  for  high  elevation  angles. 

In  either  case  the  time  available  for  CFLOS  is  higher  for  lower  closure  rates  and 
therefore  the  probability  of  obtaining  the  CFLOS  is  higher  (i.e.,  there  is  more 
opportunity) . 

The  Cloud  Characteristics 

Simple  logic  dictates  that  there  is  a direct  relationship  between  the  probabil- 
ity of  having  a CFLOS  and  the  probability  of  having  those  cloud  conditions  which 
preclude  having  a CFLOS.  Predicting  the  probability  of  having  a CFLOS,  therefore, 
requires  having  "sufficiently  detailed  cloud  information."  This  information  must 
be  merged  with  the  other  factors,  e.g,,  dynamics  and  boundary  conditions,  involved 
in  the  problem. 

For  the  problem  being  discussed,  "cloud-base"  Is  part  of  the  required  informa- 
tion. This  report  will  show  that  information  regarding  the  amount  of  sky  covered 
at  a particular  altitude  and  the  3-dimensional  geometry  and  spatial  distribution  of 
the  individual  elements  will  also  be  required.  What  is  yet  uncertain  is  the  practi- 
cal meaning  of  "sufficiently  detailed"  information  and  if,  and  how,  that  meaning 
changes  as  the  problem  variables  change. 

Historical  records  of  surface-based  weather  observations  can  be  used  to  obtain 
information  regarding  the  probability  that  a particular  sky  coverage  will  be 
experienced  at  a particular  altitude;  the  detail  available,  however,  varies  inversely 
with  th*.:  altitude.  Satellite  data  may  eventually  ue  useful  in  adding  more  detail 
at  higher  altitudes,  but  currently  the  information  is  limited  to  the  lowest  8000 
feet.  Little  or  no  information,  and  no  information,  and  no  historical  records  exist 
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relative  to  the  "individual  element"  characteristics.  Therefore,  one  must  model  the 
individual  element  characteristics  and  use  historical  data  bases  to  obtain  the  sky 
cover. 

The  authors  have  yet  to  determine  how  best  to  utilize  historical  data  since  the 
sensitivity  of  answers  to  relatively  unknown  cloud-field  characteristics,  such  as 
the  distribution  and  3-D  geometry  of  the  individual  elements,  are  unknown.  How  will 
the  Interrelationship  of  cloud  base  and  cloud  cover  as  they  relate  to  the  dynamics 
of  the  problem,  affect  the  answers?  These  ore  specific  areas  in  which  further 
analysis  is  required. 


SECTION  C — A DYNAMIC  CFLOS  COMPUTER  PROGRAM 


The  computer  program  presented  in  Appendix  D models  a situation  in  which  two 
points  are  separated  in  space  in  the  presence  of  a cloud  field  to  help  study  dynamic 
CPLOS  problems.  The  lower  point  13  begins  moving  from  its  initial  position  at  a 
fixed  heading  and  speed.  The  higher  point  A is  introduced  at  a given  vertical 
separation  and  slant  range  and  a preselected  azimuth  to  the  right  of  B's  path. 

Point  A follows  some  path,  usually  the  direct  intercept  course,  and  the  line-of- 
sight  between  A and  13  is  examined  each  0.2  seconds  of  simulated  time  until  a pre- 
viously determined  slant  range  is  reached.  B is  then  reset  to  its  original  position 
and  the  some  process  is  repeated  except  that  A approaches  from  the  preselected 
azimuth  to  the  left  of  B's  path.  B is  then  reset  to  Its  original  position,  a new 
heading  is  established  30°  to  the  left  of  the  original  heading  and  the  process 
starts  anew.  The  scenario  Is  repeated  until  a total  of  12  headings  (360°)  have  been 
examined.  The  "statistics"  are  then  summarized  for  the  2*1  approaches  involved. 

Some  Limltatlons/Wcaknesscs  of  the  Program: 

a.  CLDNO,  one  of  the  two  cloud  models,  represents  a mid-continental,  u.iiform- 
ourface  area,  carly-aftornoon,  summertime  situation.  This  situation  is  relatively 
simple  to  model.  CLDNO  produces  uniform  circular  elements,  the  centers  of  which  ore 
placed  randomly  in  the  X,Y  plane.  Coverage  is  determined  by  the  X,Y  projection  of 
the  individual  elements  onto  the  plane.  The  base  height  and  thickness  can  be  set 
but  are  uniform  for  all  elements. 


The  random  placement  of  the  elements  limits  the  total  coverage  in  the  applica- 
tion of  the  model.  Theoretically,  one  could,  by  careful  packing,  achieve  a coverage 
of  about  78#  (see  Appendix  C).  In  running  the  program,  a practical  limit  is  reached 
at  about  *15#.  The  program  CLDOVR,  developed  because  of  this  practical  limit,  allows 
for  some  restricted  overlap  of  individual  elements  and  seems  to  have  a practical 
limit  in  coverage  of  about  70%.  With  larger  fractional  coverages  the  physical  model 
more  nearly  represents  a late  afternoon  situation  but,  because  of  the  common  tops 
used,  is  less  than  an  adequate  model.  Modification  of  the  model  will  cont.nue  to 
overcome  the  "coverage,"  "uniform  top,"  and  "single  layer"  limitations/weaknesses. 

b.  The  coverage  produced  in  not  "sky-coverage"  but  rather  "earth-coverage." 

The  former  requires  a projection  of  the  cloud  elements  against 
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while  the  latter  is  obtained  by  projecting  the  cloud  bases  against  the  assumed  flat, 
underlying  surface. 

c.  The  chase  models  are  limited  in  their  applicability.  For  example,  oie  would 
not  use  either  model  to  simulate  a "dog  fight."  Both  models  currently  require  that 
A be  moving  faster  than  B. 

d.  The  assessment  of  the  line-of-sight  takes  place  every  0.2  seconds  >f  simu- 
lated time.  This  time-step  may  be  inappropriate  for  some  applications  (se»  Appen- 
dix B). 

Questions  Addressed  by  the  Program 

Given  a particular  set  of  conditions,  the  program  is  designed  to  answer  three 
questions: 

Question  1.  What  is  the  probability  of  having  a CFLOS  "at  any  instant?" 

The  line-of-sight  between  A and  B is  examined  at  each  A,B  position 
until  A and  B approach  within  S of  each  other.  The  number  of  times 
the  line-of-sight  is  free  of  clouds  (successes)  is  divided  by  the 
number  of  times  it  is  examined  (attempts).  The  results  are  accumulated 
for  each  of  the  12  headings  followed  by  B and  printed  out  as  the 
"static  case"  shown  in  Figure  5» 

Question  2.  What  is  the  probability  of  having  a CFLOS  which  lasts  for  at  least 
t -units  of  time  at  any  time? 

A time,  t,  is  chosen  as  an  integral  multiple  of  0.2  seconds.  That  is, 
t=0.2N,  where  N-1,2,3, . . .M.*  The  first  N-consecutive  lines-of-sight 
are  then  examined  for  a continuous  CFLOS.  A record  is  made  of  the  re- 
sults of  that  examination.  The  first  line-of-sight  is  then  ignored 

and  the  next  N-consecutive  lines-of-sight  are  examined,  etc.  The  number 
of  times  that  N-consecutive  lines-of-sight  are  cloud  free  divided  by 
the  number  of  times  N-consecutive  lines-of-sight  are  examined  before  S 
is  reached,  forms  the  basis  for  the  answer  to  this  question.  The 
process  is  repeated  for  all  B-headings  and  the  collected  statistics 
show  up  as  "Probability  of  a CFLOS  for  the  time  indicated,"  with  t, 
in  seconds,  following  ".GE."  in  Figure  5* 

Question  3.  What  is  the  minimum  frequency  of  occurrence  of  a CFLOS  lasting  for  at 
least  t-uni'ts  of  time? 

For  each  initial  placement  of  a,B,  a binary  record  is  made  of  whether 
or  not  a CFLOS  existed  for  t-units  of  time  at  least  once.  For  any  ini- 
tial set  of  conditions  and  for  a given  cloud  situation  the  maximum 

number  of  attempts  would  be  24;  i.e.,  A approaches  B from  left  and  right 
aspects  for  each  of  12  different  headings  used  by  B.  During  each  chase, 
a CFLOS  either  did  or  did  not  exist  at  least  once  for  t-units  of  time. 


* For  the  limit  of  M,  see  Appendix  B. 
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The  output  of  this  analysis  is  found  in  the  Minimum  frequency..." 
section  of  Figure  5. 


■Wf*  1 III— ■!  *rnm  ■ m i..  .mi 

PROBABILITY  OF  A' CFLOc  - STATIC  CASE-  .2*7 

PPOBARILITv  OF  A CFLOc  fop  at  | FAST  THF  TIMt  INOICATfO 
.ge»1  .080  .GF.?  .nn7 

MINIMUM  FMfOIIfnCv  OF  ncruPPFNCr  OF  A CFLOS  FOP  AT  LEAST  THF  TIMt  INDICATED 
OMF  SECOND  6/2 A two  SECONDS  1/2* 


A«OVF  INFO  RASFO  ON  - R ASF  = R.S.TOPS  = 11. S .RADIUS  * A.O  .COVERAGF  * 

PROBABILITY  OF  A CFLOs  - STATTr  CASE-  .2?« 

PRObARILTTv  OF  A CFLOc  FOR  AT  1 FAST  THF  TIMt  INDICATfO 
.OF. 1 .*51  .OF.?  o.ooO 

MINIMUM  FRfOUfmCy  OF  nCCUPoFNCr  OF  A CFLOS  FOR  at  I East  THF  TIMt  INDICATED 
OK'F  S>FrONO  5/24  TWO  SECONDS  0/24 


AROVF  INFO  BASFO  ON  - BASF  = S.S.TOPS  = 14.5  .RADIUS  s 4.0  .COVERAGE  * 

PROBABILITY  OF  A CFLOc;  - STAVir  CA^F"  ,356 

PROBABILITY  OF  A CFLOc  for  at  | cast  THF  TIMt  INDICATfD 
.OF.  1 .121  .OF.?  .0<>7 

MINIMUM  FRfOUfnCY  OF  oCCURRF'OCr  OF  A CFLOS  FOR  AT  LEAST  THF  TIME  TNOTC^tfO 
OMF  bFCOND  10/24  TWO  SECONns  1/24 


AR0VE  INFO  RASFO  ON  - RASE  * 7.0*TAf'1;  * 14.5  .RADIUS  s 4.0  .COVERAGE  * 

PROBABILITY  of  A CFLOc  - STATIC  CASE-  .356 
PROBABILITY  OF  A CFLOc  FOR  AT  I FAST  THF  TIMt  INDICATfO 
, OF  * 1 .121  .OF.?  .« 07 

MINIMUM  ERfOUCMCy  OF  orCURRFNCF  OF  A CFLOS  FOR  AT  LEAST  THF  TIMt  INOICATFD 
ONE  SECOND  10/24  Two  SECONDS  1/24 


ABOVE  INFO  BASFO  ON  - RASE  = y.O.TOPS  = 18.3  .RADIUS  * 4.0  .COVERAGE  * 

PBOHARILITy  of  A CFI.Oc  - STATIC  CASE-  .310 

PROBABILITY  OF  A CF.Os  FOR  At  i FAST  THF  TIMt  INDICATfO 
.OF. I .174  .OF..?  .097 

M !i ‘ 1 MUM  FRfQUFNCY  OF  pCCURPFNC-'  OF  A CFLOS  FOR  AT  t EAST  THF  TIMt  INDICATFO 
ONF  bF<  OND  9/24  TWO  SECONDS  3/24 


AROVF  INFO  BASF D ON  - RASE  = 4.0. TOPS  = 9.0  .RADIUS  = 1.0  .COVERAGE  « 

PROBABILITY  OF  A CFLOc  - ST  AT  I*  CASF-  »27fl 

PROBABILITY  OF  A CFLOs  FOR  AT  | FAST  THF  TIMt  INDICATfO 
,GE» 1 .136  .GF.?  .069 

MINIMUM  FRFQUFNCY  OF  pCCUPPFMC.f  OF  A CFLOS  FOR  AT  LEAST  THF  TIMt  INDICATED 
ONF  SECOND  0/24  TWO  SECONDS  2/24 


AHOVF  fNFO  BASFO  ON  - RASE  = 4.0. TOPS  = 11.5  .RADIUS  = 1.0  .COVERAGF  * 


.435 


.435 


.435 


.435 


.245 


Figure  5-  Sample  Output  of  the  Program. 
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Some  Model  Results 

Figures  ha  and  6b  show  scaled  drawings  of  one  quadrant  of  a cloud  field  gener- 
ated by  CLOUDNO.  The  light  dashed  lines  represent  B's  path,  which  is  actually  at  c 
height  which  is  below  the  cloud  base.  The  heavy  dashed  lines  are  the  projections  of 
A 1 s path  into  the  plane  of  B.  The  heavy  solid  lines  ore  the  actual  paths  followed 
by  A.  See  Appendix  D for  a brief  description  of  the  HUNT  and  HUNTD  chase  models, 
in  effect  then.  Figures  6a  and  6b  are  pictures  of  a "run."  The  lines-of-sight 
which  arc  examined  arc  those  lines  which  would  connect  the  positions  of  A and  B at 
any  time. 


Figure  6a.  A Visual  Kepresentation  of  Figure  6b.  A Visual  Representation  of 

CLOUDNO  and  HUNTD.  CLOUDNO  and  HUNT. 

Table  1 presents  the  statistics  associated  with  many  such  runs.  The  three  sets 
of  statistics  presented  in  the  table  are  related  to  those  questions  discussed  in  the 
previous  section.  For  a given  run  there  are  approximately  600  "instantaneous- 
attempts"  and  *150-550  one-second  "dynamic-attempts."  The  "minimum  frequency"  values 
cannot  exceed  2*1  since  there  are  only  2*i  opportunities  per  run  for  having  a CFLOS 
for  the  time  span  indicated  at  least  once.  There  is,  of  course,  only  one  value  of 
"earth-cover"  for  a given  run. 

The  two  sets  of  statistics  associated  with  A and  E were  derived  from  the  same  9.2 
nonoverlapping  cloud  distributions.  The  distributions  were  established  based  on  a 
*!00-mcter  cloud  radius.  Since  the  cloud  elements  do  not  overlap,  one  could  simply 
change  the  radius  and  calculate  the  new  coverage.  Most  of  the  data  associated  with 
the  B and  F statistics  are  derived  from  these  same  distributions  by  simply  changing 
the  thickness  of  the  cloud  elements.  About  one-half  of  the  C and  P statistics  use 
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these  distributions,  so  there  are,  at  most,  30  different  cloud  element  distributions 
certainly  not  a representative  sample  of  the  infinite  variety  of  distributions 
which  could  be  produced  either  in  the  real  world  or  in  the  model.  Nonetheless, 

Table  1 can  be  used  to  gain  insight  into  dynamic  CFLOS  problems. 

The  content  of  Table  1 is  associated  with  a very  specific  set  of  parameters ; 
i.e.,  specific  values  for  V^,  2^,  S,  and  a.  The  authors  wish  to  stress  the 

point:  Dynamic  CFLOS  problems  tend  to  be  complex  and  unique  and  it  will  be  diffi- 
cult to  generalize  solutions  for  large  ranges  of  the  associated  parameters. 

The  "Earth-Cover"  for  Cases  A,  B,  C,  and  D are  essentially  the  same.  Examina- 
tion of  the  associated  "instantaneous"  and  "dynamic"  probabilities  indicates  that 
the  values  may  be  more  sensitive  to  changes  in  cloud-base  height  than  to  changes  in 
cloud  thickness.  This  sensitivity  relationship  is  even  more  evident  in  the  "Mini- 
mum Thickness..."  data.  Comparing  E and  F to  C leads  to  the  tentative  conclusion 
that,  at  least  for  some  combinations  of  decreased-coverage/increased-cloud-base, 
the  latter  will  be  the  more  significant.  Finally,  the  probability  of  having  a 
finite -time  CFLOS  is  less  than  the  probability  of  having  a CFLOS  instantaneously. 
This  follows  from  the  fact  that  a low  instantaneous  value  indicates  a frequently 
obscured  line -of- sight. 

Experience  indicates  that  the  distribution  of  the  individual  elements  is  a sig- 
nificant variable  and  that  its  significance  increases  as  earth-coverage  decreases. 
This  seems  to  be  particularly  true  for  earth-coverage  less  than  about  *10#.  The 
authors  also  tend  to  believe  that,  for  a given  coverage,  "cloud  radius"  is  not  a 
significant  variable;  although  the  elements  are  smaller,  they  are  more  numerous  and 
the  increased  population  effectively  cancels  the  effect  of  the  decreased  size. 

Changes  in  V^/V^  are  significant,  as  can  be  seen  by  a comparison  of  Table  1 and 
Table  2.  The  same  specific  parameters  are  used  for  both  tables  except  that,  in 
Table  2,  V^/V^  ~ 3,  and  10  completely  different  distributions  were  used.  Case  1 of 
Table  2 and  Case  A of  Table  1 differ  only  slightly  in  mean  earth-coverage  and  may  be 
compared  directly.  Being  cautious  to  note  the  differences  in  cloud  characteristics, 
one  may  also  compare  Case  2 to  Cases  C and  D,  and  Case  3 to  Cases  E and  F.  The 
generally  higher  values  of  Table  2 would  seem  to  stem  primarily  from  the  lower  speed 
ratios;  but  one  must  bo  cautious  of  such  a generalization. 


Table 

CFLOS 

Probabilities  for 

Various  Modeled  Cloud  Situations  when 

VVB 

= 3, 

VB 

2*30  m/soc 

, S = 001 

meters,  /.^  = 

100  meters. 

and  e 22°. 

CFLOS  Results 

Cloud  Data 

Pr 
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Tabic  3 presents  selected  calculations  of  Z ^ which  can  be  used  to  help  interpret 
the  statistics  of  Tables  1 and  2.  (Linear  interpolation  of  the  values  will  be 
accurate  to  within  about  l£.)  For  the  conditions  in  Table  1,  Zc  548  meters. 

Since  Zg  for  Tables  1 and  2 was  150  meters,  any  cloud  base  above  698  meters  for 
Table  2,  or  above  795  meters  for  Table  1,  would  assure  a CFLOS  for  one-second  at 
least  once.  Therefore,  the  minimum  frequency  values  of  Case  C and  Case  2 are  essen- 
tially predictable  a priori.  Note  also  that  if  A approaches  B from  behind  with  a 
relative  azimuth  of  165°  VA/Vy  “ 3,  0 * 20°,  the  minimum  frequency  column  of  Case  D 
could  be  "predicted"  (i.e.,  calculated)  a priori. 

Table  3*  Z^,  Values  for  Selected  Conditions 

Vg  = 250  m/sec,  t = 1 sec,  S *»  501  meters. 

Speed  Ratio:  VaAr  * 3 Speed  Ratio:  VA/Vg  = ^ 

Initial  Relative  Azimuth 
Elbv.  ■ — — 


Angle 

15° 

30° 

45° 

165° 

150° 

135° 

15° 

30° 

450 

165° 

150° 

135' 

10° 

258 

252 

244 

175 

178 

183 

301 

296 

288 

219 

222 

227 

20° 

503 

492 

476 

347 

353 

363 

589 

579 

564 

434 

440 

450 

CO 

O 

0 

723 

710 

688 

514 

522 

535 

850 

837 

817 

641 

649 

664 

SECTION  D — SUMMARY  AND  CONCLUSIONS 

Current  techniques  for  predicting  the  probability  of  a CFLOS  are  designed  for 
problems  which  are  instantaneous  in  nature  and  involve  stationary  points.  Dynamic 
CFLOS  problems  Involve  either  the  movement  of  one,  or  both,  of  the  points  b tween 
which  the  line-of-sight  is  to  be  assessed  or  a time  during  whicn  the  line-of-sight 
is  to  be  assessed,  or  both  movement  and  time  may  be  involved.  Hence,  current  CFLOS 
prediction  techniques  can  not  be  used  for  dynamic  CFLOS  problems.  To  assess  the 
potential  usability  of  those  airborne  electro-optical  systems  which  are  effectively 
blocked  by  clouds,  new  techniques  must  be  developed.  "Modeling"  is  the  approach 
which  the  authors  have  used  as  a preliminary  step  in  that  direction. 

The  computer  program  discussed  herein  indicates  that  dynamic  CFLOS  problems  tend 
to  be  complex,  unique,  and  not  amenable  to  generalization.  To  solve  a dynamic  CFLOS 
problem  one  requires  a detailed  knowledge  of  the  initial  positions  and  subsequent 
velocities  of  the  two  points  between  which  the  line-of-sight  is  to  be  assessed, 
detailed  knowledge  of  the  boundary  conditions  related  to  the  space/time  dimensions  of 
the  problem,  and  detailed  knowledge  of  the  cloud  field  characteristics,  including 
the  3-dimenalonal  geometry  and  spatial  distribution  of  the  individual  cloud  elements. 
The  dynamics  of  a problem  might  be  definable  through  a knowledge  of  the  operational 
characteristics  of  the  aircraft  involved.  The  boundary  conditions  and  time  incre- 
ments to  be  used  in  the  line-of-sight  assessment  might  stem  from  the  characteristics 
of  the  electro-optical  system.  These  are,  essentially,  nonmeteorological  problems. 
The  meteorologist  will  have  to  define  the  cloud  field  characteristics. 
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Practically  speaking,  extensive  knowledge  of  the  required  cloud  field  character- 
istics is  limited  to  "cloud  amount  vs  altitude"  and  even  that  is  limited,  essentially, 
to  the  lowest  8000  feet.  Satellite  data  may  be  useful  in  extending  our  knowledge  of 
this  parameter  at  higher  altitudes.  Detailed  information  relative  to  the  geometry 
and  spatial  distribution  of  the  individual  elements  is  extremely  limited  to  nonexist- 
ent. To  determine  what  degree  of  resolution  will  be  required  and  how  much  intelli- 
gence can  be  derived  from  historical  records  will  require  further  insight  into  the 
nature  of  dynamic  CFLOS  problems  and  further  analysis  cf  the  data  base. 

The  computer  program  presented  herein  can  be  a useful  tool  in  developing  insight 
into  the  variables  associated  with  dynamic  CFLOS  problems.  In  particular,  it  can  be 
used  to  examine  the  sensitivity  of  the  answers  to  changes  in  the  input  parameters. 
Being  careful  to  observe  the  artificiality  of  the  model,  it  can  also  be  used  to 
obtain  what  might  best  be  termed  "bpll-park"  estimates  of  real  world  situations. 

For  example,  for  situations  similar  to  those  discussed  in  this  report,  it  would  seem 
that  the  probability  of  having  a cloud-free  line-of-sight,  either  instantaneously  or 
for  as  short  a period  as  one  second,  will  be  low  even  in  fair-weather  conditions. 
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Appendix  A 

DIRECT  INTERCEPT  MODEL 


Problem:  Given  the  initial  positions  and  speeds  of  A and  B and  a level  route  of 
flight  at  known  azimuth  for  B,  find  the  distance  A must  travel  to  inter 
cept  B. 


AQ,  Bc  * Initial  positions  of  A and  B 

o = Elevation  angle  from  B„  to  A 

o o 

= Azimuth  of  the  path  to  be 
L>  traversed  by  B 

0,  * Azimuth  of  A„ 

A o 

v = Angle  between  AqBq  and  B0B' 


Figure  A-l.  Direct  Intercept  Model. 

Let  the  following  definitions  apply: 

(XA,  Ya,  2a),  (0,  0,  Za):  x,  y,  z coordinates  of  AQ  and  Bq 

C : Distance  W 

D : Distance  TIT 
o o 

Da,  Dq  : Distance  to  be  traveled  by  A and  B 
VA,  VQ  : Constant  speeds  of  A and  B 


Begin  by  solving  for  Cos  y.  From  Figui*e  A-l  and  the  definitions  given, 

YA  = D Cos  a Sin  $A  (A-l) 

XA  = D Cos  o Cos  &A  (A-2) 

W s h Sec  °B 


(A-3) 
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From  the  law  of  cosines  we  have: 

C2  - D2  + (E^)2  - 2D  (SJF)  Cos  Y 
and,  substituting  from  Equation  (A— 3 ) 

C2  - D2  + X2  Sec2  eB  - 2DXa  Sec  f>B  Cos  y (A -4) 

But,  note  that 

C2  » (ZA  - ZD)2  + (Ya  - XA  Tan  6B)2  (A-5) 

Hence,  one  can  set  Equation  (A-h)  equal  to  Equation  (A-5)>  carry  out  the  multipli- 
cation indicated  in  the  second  term  of  the  right  side  of  Equation  (A-3),  and 
rearrange  and  combine  terms  to  obtain 


Cos  y 


XA  f YA  Ta»  °B 
D Sec  SB 


(A-6) 


Substituting  Equations  (A-l)  and  (A-2)  for  YA  and  XA,  respectively,  and  rearranging 
and  simplifying  terms  results  in  the  equation 


Cos  Y = Cos  Q 

(Cos 

0A  + Sin  SA  Tan  Ag)  Cos  °B 

(A-7) 

which  will  be  generalized  for  use 

as 

Co 

s Y * Pj_ 

(A-8) 

By  definition,  DA  = and  Dg  = 

VBt. 

Therefore, 

db 

-d  !s 
*'« 

(A-9) 

From  the  law  of  cosines  we  have: 

D2  » D2  + - 2DD0Cos  Y (A-10) 

Substituting  from  Equations  (A-8)  and  (A-9)  and  rearranging  terms 

da  - ;§)  + rr  pi  - 1,2  - 0 <‘-n) 

VA  A 

from  which,  using  the  general  solution  to  quadratic  equations,  and  requiring  that 


D 


± 


+ 1 


(A -12) 


1 
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Or,  for  VD  = VA 


DA  = 


(A-13) 


and,  physically,  only  the  positive  radical  of  Equation  (A-12)  applies. 

In  the  text  the  problem  was  stated  assuming  \fA  > Vg.  The  computer  program 
currently  requires  that  VA  be  greater  than  Vg.  However,  Equations  (A-12)  and  (A-13) 
clearly  indicate  that  DA  has  solutions  which  are  not  restricted  by  VA  > Vg.  Three 
possibilities  will  be  examined:  VA  = Vg,  > Vr<>  and  < Vg. 


Case  1:  VA  = Vg 

From  Equation  (A-8),  the  range  of  is:  - 1 * F^  < 1.  Physical  reasoning 

dictates  that  F-^  < 0 be  eliminated  since  A cannot  overtake  B when  approaching  from 
behind  B's  initial  position.  Therefore,  only  0 * F^  * 1 is  considered. 

From  'iquation  (A-13)  and  the  definition  of  F-^ 

§-  = 2 Cos  y (A-1^4) 

!)A 


from  which  solutions  may  be  obtained  when  VA  = Vg  for  all  choices  of  y such  that 
< Y < J . 

V 

Now  examining  rr^  = K where  K i*  1 and  using  Equation  (A-12),  we  obtain 
VA 


da=  D 


- KFl  ± / K2  (F2  - 1)  + 1 


1 - Kc 


(A-15) 


Given  that  K ^ 1,  Equation  (A-15)  only  requires  that 

K2  (F2  - 1)  * - 1 (A-16) 

2 2 

Clearly  any  negative  contribution  must  come  from  F-^  - 1.  This  dictates  that  K must 
be  restricted  to  certain  maximum  values  to  insure  that  the  condition  of  Equation 
(A-l6)  is  not  violated.  Table  A-l  and  Figure  A-2  illustrate  the  values  as  a function 
of  y for  the  range  0 s y < 90°. 


Tables  A-l 


kmax 

00 

5.75 

2.9: 

Y 

0 

10 

20 

Maximum  K-Value 

2.00  1.5b  1.31 
30  40  50 


1.15  1.0(5  1.02 

60  70  80 
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Case  2:  V.  > Vr, 

For  this  case  K < 1 and  therefore 
the  condition  imposed  by  Equation  (A-l6) 
is  always  met. 

Case  3:  VA  < V0 

For  this  case  K > 1 and  therefore 
K must  be  restricted  to  the  range 
1 < K < K 

max 

Finally, 

da  = df2  (A-17) 

where  Fg  is  defined  according  to  the 
following: 

Case  1:  Fg  «*  when  VA  = V0 


Y 


Figure  A-2.  Violation  Area  for 
K and  Y. 


Case  2: 


when  V.  > V., 
A B 


Case  3: 


F2  = 


KF1  + J K2  (F2  - 1)  + 1 


subject  to  both  of  the  following  conditions: 

a.  K2  (F2  - 1)  * - 1 

b.  K > 1 


Note  that  Case  2 is  the  case  discussed  in  the  text. 
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Appendix  B 

ANALYSIS  OF  THE  CRITICAL  RELATIVE  CLOUD  BASE 

Problem:  Determine  the  height,  Zc,  above  the  plane  of  B,  above  which  the  existence 
of  clouds  will  have  no  effect  on  the  probability  of  having  a CPLOS  between 
A and  B for  At  units  of  time  at  least  once  before  A and  B approach  within 
distance  S of  each  other,  given  that  A pursues  B in  accordance  with  the 
direct  intercept  model. 


0 : Azimuthal  separation  of 

A^,  and  B^  (relative 
o o ' 

azimuth) 

0 : Elevation  angle  from 

K to  A 
o o 

D : Distance  A B 
o o 

Z : Base  height  of  the  low- 
est layer  of  clouds 

P : Intersection  of  the 
paths  of  A and  B 


Figure  B-l . Critical  Cloud  Base. 


Let  the  following  definitions  apply: 


Brv»  D>  Da>  db*  zaj  ZB: 


C" 


As  defined  in  Appendix  A 
Angle  between  A0P  and 

Distance  from  A„  to  Z along  the  path  of  A 

Distance  from  the  intersection  of  A's  path  with  the  plane 
of  Z to  that  point  where  A and  B come  within  distance  S 
of  each  other 
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A B 

Lg,  Lg:  Distances  which  have  yet  to  be  traveled  by  A and  B from 
that-  point  where  A and  B are  within  distance  S of  each 
other  to  the  point  P. 


Noting  from  the-  definitions  that: 

DA=  LC 

and  therefore 


+ 4 + L 

- Lc  - 4 


and 


Cos 


(B-l) 


from  which 


From  similar  triangles 


By  definition,  At  = ~ . 
obtain:  A 


Combining  this  with  Equations 


(B-2) 

(B-3) 

(B-l),  (B-2),  and  (B-3),  we 


from  which 


D' 


+ <ZA 


AtV 


A\ 

/ 


(B-4) 


Hence  the  cloud  base  relative  to  the  height  of  B is 


AtV 


AN 


ZC  - Z - h - <ZA  - ZB>  (!  + -D -J 
Substituting  from  Equation  (A-17),  Appendix  A,  for  DA  yields 

AtV. 


ZC  - <ZA  • ZB>  (l  + IT*-) 


and,  since  (ZA  - Zg)  = D Sin  a 


Zc  = Sin 


/ *tVA' 


(B-5) 


(B-6) 


Hence,  given  VA  and  the  definition  of  Fg  (see  Appendix  A),  one  can  choose  a 
slant  range  (S)  aid  At  of  interest  and  calculate  the  height  above  which  cloud  bases 
will  have  no  impact  on  the  probability  of  having  a CFLOS  between  A and  B at  least 
once  before  A and  B approach  within  S of  each  other. 
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Clearly,  At  will  not  be  meaningful  unless  it  is  chosen  such  that 


and  therefore 


At  < 


DF„ 


At  < 


In  the  analysis  procedure  described  in  the  text  At  ? treated  as 


(B-7) 


At  » gN 

where  g is  the  incremental  time-step  used  in  the  analysis  of  the  line-of-sight  and 
N has  a maximum  value  of  M such  that  the  total  time  involved  in  the  chase  is  not 
exceeded.  Clearly,  g only  makes  sense  when  it  is  some  fraction  of  At,  but  what 
fraction  should  it  be?  A value  of  0.2  seconds  was  chosen  for  the  program  based  on 
a subjective  analysis  of  the  interrelationship  of  the  cloud  characteristics  (radius 
and  thickness)  and  the  speeds  typically  involved.  The  authors  have  not  yet  examined 
the  sensitivity  of  th*'  results  to  changes  in  "g-values"  but  this,  too,  is  an  area 
which  needs  study.  In  any  case,  one  should  not  choose  g completely  arbitrarily. 
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Appendix  C 

CLOUD-COVER  LIMITATION  OP  A UNIFORM  DISTRIBUTION  MODEL 


Problem:  Determine  the  maximum  fractional  coverage  obtainable  from  placing  uniform 
circular  elements  into  the  large  circular  area  such  that  each  element  is 
totally  contained  within  the  larger  area  and  no  elements  overlap. 


r = Radius  of  cloud  elements 


R = Radius  of  entire  area 

R,  = Radius  of  ith  locus  of  cloud 
1 centers  where 


Ri  < Ri+1 


Figure  C-l.  A Uniform  Distribution 
Model. 

Lot  the  following  definitions  apply: 

L^  = Circumference  of  i^  locus 
N^  a Number  of  elements  on  L^ 

N c Total  number  of  elements 
imnx  = Maximum  number  of  loci 

Case  1:  Condition  0 <r  < < R 

The  circumference  of  the  l^1  locus  is 

I.l  = 2ttR1  (C-l) 

The  radius  of  the  i^11  locus  is 


Ri  = 2ri 


(C-2) 


therefore  L^  *»  l(nri 

If  the  number  of  elements  on  a given  locus  must  be  an  integral  value,  then 


(C-3) 


Nt  « 2ni 
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The  maximum  number  of  loci  which  can  be  allowed  for  a given  cloud  radius  will  be 


R 

27- 


Prom  Equation  (C-2),  we  note  that  a hole  of  radius  r has  been  left  at  the  center 

1) 

of  the  area.  Hence,  if  =>  K,  and  K is  truncated  to  an  integer,  then  i^^.  «=  k-1 
and 

K-1 

N < 1 •+•  ^ 
i.Hl 

and  from  N.  *■  2ni 

1 K-1 

N * 1 + Sn  \ i (C-'i) 

J=1 


and  testing  various  values  of  K,  the  authors  found  that 


N * 1 + n K(K-l) 

The  fraction  of  the  area  covered  by  clouds  is 

p = Nm£  = Nr! 
TiR**  R2 


and  therefore 


Finally,  we  note 


4kc 


lim  P ■ £ - . 78!54 

k-*» 


(C-r>) 


In  practice,  the  authors  have  observed  that  randomly  placed  elements,  which  are  not 
allowed  to  overlap,  rarely  exceed  5 6#  coverage. 

Case  2;  0 < r < R 

The  solution  of  r **  P and  N = 1 with  the  element  centered  at  the  center  of  the 
circle  is  a trivial  but  valid  solution.  Another  solution  is  R for  which  P » jp 
Hence,  one  can  obtain  fractional  coverage  in  excess  of  ^ for  a single  element 
centered  at  the  center  of  the  circle  whose  radius  meets  the  condition 

»/?<«•<  R 
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Appendix  D 

LISTING  OF  A DYNAMIC  CKLOS  COMPUTER  PROGRAM 


The  computer  program  discussed  herein  can  be  thought  of  in  terms  of  two  separate, 
but  relatable,  parts:  Driver  programs  and  subroutines.  Driver  programs  are  used  to 
define  the  values  which  the  variables  are  to  have  in  the  subroutines  and  to  model 
the  "evolution"  of  the  dynamic  encounters.  The  currently  modeled  evolution  for  all 
driver  programs  is  as  described  in  Section  C of  this  report.  Alternative  models 
might  call  for  a single  azimuthal  approach  (e.g.,  15°  to  the  left  of  D,  as  a single 
choice),  but  repeated  with  higher  resolution  (e.g.,  reanalysis  for  incremental 
changes  of  10°  in  B's  heading).  A better  understanding  of  the  overall  program  can 
be  obtained  by  operating  one  of  the  driver  programs  described  below. 

Driver  Programs 

I3CFL0S  A basic  driver  program  which  utilizes  a nonovc slapping  cloud  field,  a direct 
intercept  chase  model  and  azimuthal  variations  of  ± 25°  from  B's  heading. 
Summaries  are  made  of  the  line-of-sight  statistics  collected  for  the  12 
headings  taken  by  B.  BCFLOS  could  be  run  with  an  overlapping  cloud-field. 
(Table  D-l.) 

IDS0AZ  A somewhat  more  involved  program  which  is  to  be  used  only  with  a nonover- 
lapping cloud-field.  Reanalysis  of  the  line-of-sight  is  performed  as 
changes  are  introduced  in  cloud  base,  thickness,  and  radius  values,  for  a 
given  cloud-element  distribution.  (Table  D-2.) 

L0S1AZ  A program  designed  specifically  for  the  overlapping  cloud- field.  Users 
will  note  that  the  only  significant  difference  between  LOS0AZ  and  L0S1AZ 
is  that  in  the  latter,  the  cloud-radius  is  not  varied.  (Table  D-3.) 

Subroutines 

Comment  cards  have  been  inserted  liberally  in  all  subroutines.  However,  a brief, 
functional  description  of  each  subroutine  appears  below.  An  examination  of  the 
listing  in  Table  D-^l  will  reveal  that  many  "PRINT"  statements  have  been  rendered 
inactive  by  a "C"  in  Column  1.  These  statements  are  rarely  used  after  familiarity 
is  gained  with  the  overall  program.  The  inactive  statements  have  been  retained  in 
the  listing  to  provide  potential  users  with  tools  which  can  be  used  to  follow  the 
detailed  flow  of  the  CFL05  analysis.  Reactivation  of  all  "PRINT"  statements  will 
produce  a computer  printout  about  *1  inches  thick.  If  wholesale  reactivation  is 
desired,  we  recommend  using  BCFLOS  as  the  driver  program  since  it  is  set  up  to  pro- 
duce only  a single  cloud  distribution  and  to  "fly  through"  only  once. 

SET  Initializes  required  values  within  other  subroutines. 

CLDNO  Contains  only  otic  entry:  Entry  CLOUDNO.  CLOUDNO  establishes  the  positions 
of  the  centers  of  each  cloud  and  a percentage  coverage.  The  cloud  center 
coordinates  arc  generated  using  a random  number  generator.  Cloud  elements 
will  not  overlap. 
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CLDOVR 


HUNTF 


FREEF 


MAPF  & 
CATCH 


Contains?  a single  entry,  CLOUDNO,  which  allows  a newly  introduced  element 
to  overlap  one  existing  element  up  to  a limit  equal  to  the  cloud  radius. 
This  entry  should  be  used  for  coverages  in  excess  of  about  45#. 

Contains  two  different  entries: 

a.  Entry  HUNT  produces  the  positions  of  two  moving  objects  (A,B)  whose 
initial  positions  have  been  specified  in  the  driver  program.  The  dynamics 
of  this  entry  are  such  as  might  correspond  to  an  IR  Sensor  tracking  a 
moving  object.  In  essence,  A continually  heads  toward  the  current  position 
of  B. 

b.  Entry  HUNTD  is  a slightly  more  sophisticated  chase  routine  (see 
Appendix  A)  in  which  the  point  toward  which  A must  move  to  intercept  B is 
predetermined,  based  on  the  assumption  that  B will  continue  on  a straight- 
line  course  at  a constant  speed. 

Output  from  both  of  these  entries  includes  x,y,z  coordinates,  the  slant 
range  distance  between  A and  B,  and  the  elevation  angle  from  B to  A,  all 
as  a function  of  time.  The  parameters  are  no  longer  calculated  after  a 
predetermined  time  and/or  slant  range  has  been  reached.  LIMA  and  LIMB  are 
used  to  indicate  that  the  elevation  angle  has  exceeded  certain  boundaries, 
in  the  case  of  LIMA,  or  that  both  vehicles  are  below  the  base  of  the  cloud 
(LIMB).  Although  not  currently  used,  these  indicators  are  useful  diagnos- 
tic tools. 

Contains  only  Entry  FREE.  FREE  uses  the  positional  information  derived 
from  either  HUNT  or  HUNTD  to  examine  the  line-of-sight  to  ascertain  whether 
or  not  it  intercepts  any  clouds.  If  it  intercepts  any  cloud  it  is  deter- 
mined to  be  obscured;  otherwise  a CFLOS  "hit"  is  recorded. 

These  are  the  "bookkeeping"  and  "statistics  gathering"  subroutines. 
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Tabic  D-l.  Listing  of  BCFLOS. 


PROGRAM  BCFL OS* INPUT ,OUTPUT<  A 1 

.C _IH1S„PRQGRAM„MAS,  .DES1GNED_TQ_EXAM1NE._THE™UME-.H1ST.QRY„0F.J HE.  CLOUD  ,A„  ,2 

C OBSTRUCTED  OR  CLOUD  FREE  LINE  OF  SIGHT  BETWEEN  TWO  OBJECTS  MOVING  A 3 

£ IN_THE_P.RESENC£  OF.  A.  CLOUO._R.IELO..  THE  .USER.JBUST..SREC  IF.Y..THE A_  _4. 

ikit  v i /•  Mtin  • t • rtkir  aattv  • m'  ntir  «*  run  nurr  u tun  n tir  a r 


C INITIAL  CONDITIONS  AND  OECIOf  ON  ONE  OF  TWO  CHASE  MODELS  AND  ONE  A 5 

C OF  TWO  CLOUD  MODELS  TO  BE  USED  IN  THE  ANALYSIS.  AFR1AL  COVERAGE  OF A 6. 

C THE  CLOUD  IS  DETERMINED  BY  THE  FRACTION  OF  THE  UNDERLYING  SURFACE  A 7 

X. _CQV.EREU._8Y.  .T.HE  JIRO J EC  U 0N  .QF_.T  HE_CL0U  D„B  AS  E S . A _ 8. 

COMMON  /BLK1/  XCLDX1000<, YCLDtl000<, JJ  A 9 

_COMNON_/.BLK2/_PXAX100_<tP„YAX10CKt.P.ZAX100.<»PXBXlQQ.<tR.Y.6XlQ0.<tPZ3Xl0Q  A 1.0. 

l<tTAB*100<»0A8X100<«KK  A 11 

COMMON  /BLK3/  AZA.AZB.ELtXBl  .YB1.2B1  .D1  .VA.VB A 12. 

COMMON  /BLK4/  XLNX200<,KL,L!MAX4<,L 1MBX4<  A 13 

C0MMQN_/.BLK5/_CLDB A. LA. 

COMMON  /BLK6/  J A 15 

___  _CaNN0N_/.BLK7./..BCLQ,.ICL0»PCVR. RADNOR . .A..  ,16. 

C ESTABLISH  INITIAL  CONDITIONS  A 17 

£ 2 _J  _ V ESX1CAL— SERARAU ON  Of  THE  .TWO  .OBJECTS A LR_ 

C 01  # INITIAL  SLANT  RANGE  DISTANCE  A 19 

£ EL  .§  ELEVATION  .ANGLE  FROM.  THE  SLOWER..XBC  TO  THE  EASTER  XA<  A.  20 

C OBJECT.  A 21 

£_  ..  XBlt.YU  It  ZB  I.  INITIAL  .POSITION  .OF  B.  __  A 22 

C VA»VB  • SPEEDS  OF  A AND  B.  VA  MUST  BE  GREATER  THAN  VU.  A 23 

£ RAQ_I_RAU  LUS_QF_£LQUD_tLfcKE  NLSZJLQQ* A 2A. 

C HOR  « RADIUS  OF  REGION  OVER  WHICH  CLOUDS  ARE  EMPLACED.  A 25 

£_  ._ P.CVR,.#.  PERCENT  . OF  COVERAGE  DESIRED.  A.  26 

C BCLO.TCLD  # BASE  AND  TOP  OF  CLOUD  ELEMENT S/ 100.  A 27 

. .Z#5000./3.2B,  A 28 

DI 04000.  A 29 

Y.IZ/I11 A 30. 

EL0AS!NXY<  A 31 

ELiELt.lB0./.3. 1.416  A.  32 

XBI0YBI0O.  A 33 

ZB  1 1.150.  A 34 

VAI1000.  A 35 

i ViUL25Q. A ...36.. 

HOR05OOO.  A 37 

{...  P.CVR  0.3.  A 38. 

RAD04.  A 39 

BCLO04. A ,40 

TCL0I9.  A 41 


CALL... SET.  .A.  42 

C OECIDF  ON  USING  AN  OVERLAPPING  CLOUD  FIELD  XCLOUOO<  OR  A A 43 

£ NONsOVERLAERJNG— F.LELO— XCLOUONOC.— CLOUDO — IS — RECOMMENDED— EUR A — 44. 

C COVERAGES  OF  .50  OR  GREATER.  A 45 

CALL.  CLOUDNO  A 46 

CALL  CATCH1  A 47 

, _ .CALL.  .GILO.  A.  4B. 

00  1 Jil, 12  A 49 

AZai3D.D.«J=2< — _._A 50. 

CALL  CLR  A 51 

.AZAIAZ0E25.Q  A 52 
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Table  D-l.  Listing  of  BCFLOS  (Cont'd). 


c 

OECIOE  ON  A OIRECT 

INTERCEPT  THUNTD<  OR 

TAIL  CHAS.iNG  XHUNT< 

A 

53 

c. 

^fFNARIfl. 

A 

5A 

CALL  HUNTO 

A 

55 

KKRKK-1 

A 

56 

CALL  CLRX 

A 

57 

r*l  1 FRFF 

A 

„58- 

AZA#AZB-2S.O 

A 

59 

r At  t HI  INTO 

A 

__  60- 

KKAKK-1 

A 

61 

m i riB» 

A 

A? 

CALL  FREE 

A 

63 

CALL  GIL  1 

_ A A£ 

1 CONTINUE 

A 

65 

CALL..GIL2 

A- 

—66— 

c 

***PCVR  IS  RESET  TO 

THE  AC1UAL  COVERAGE 

ATTAINED  IN  THE 

CLOUD 

A 

67 

ClIRRnilTTNF.  ^FTTING 

rnuFR  m prun  nm  a in  ^ 

THI  VAI  IIP  FflR 

PR  I NTDIIT 

A 

AH 

c 

PURPOSES.*** 

A 

69 

XQVERARCVR A.... 7.0- 

PRINT  2«  BCLD,TCLD,RAD, COVER  A 71 


2 FORMAT 

X// . 10X» 29HA80VE  INFO  BASED  ON  - BASE  f » F5. 1 .8H.TOPS  # ,F5 

A 

73 

l -1f 1 IH 

rR  Am  IK  M , F*>  . t r 1 3H  ffTU/FRAfiF  A (FA.^U" 

A 

74 

END 

A 

75- 
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Table  D-2.  Listing  of  LOS0AZ . 


PROGRAM .LOSOAZXINPUT  ,QUTPUT<  ...  - A 1 

THIS  PROGRAM  MAS  DESIGNED  TO  EXAMINE  THE  TIME  HISTORY  OF  THE  CLOUD  A 2 

.OBSTRUCT EQ-OR  CLOUD  FREE  LINE  -OF  S IGHT-.  BET  WE  EN--T  MO-OBJECTS- MOVING A—  3 

IN  THE  PRESENCE  OF  A CLOUO  FIELD.  THE  USER  MUST  SPECIFY  THE  A 

USEO  IN  THE  ANALYStS.  THIS  PROGRAM  IS  TO  BE  USED  ONLY  WITH  NON-r  A 6 

-OVERLAPPING  - CLOUD-.FI ELOS.  - ~~  A 7. 

AERIAL  COVERAGE  OF  THE  CLOUD  IS  DETERMINED  BY  THE  FRACTION  OF  THE  A 8 
-UNDERLYING-SURFACE- COVERED  BY-  THE  PROJECTION  OF.  THE.  CLOUD-BASES.-  A-  .9 
COMMON  /BLKI/  XCLDf 1000<,YCLD*1000<. JJ  A 10 

K.  TABtlOO<(DABtlOO<«KK  A 12 

-.COMMON  /BLK3/  A2A,A2B,EL,XBI  »YBI »2BI ,DI ♦ VA.VB A--13 

COMMON  /8LK4/  KLNt200<,KL»LIMAX4<,LIMBX4<  A 14 

COMMON  /BLK5/  CLDB  ..  _ _ .A..  .IS 

COMMON  /8LK6/  J A 16 

— COMMtlN-/BLiqy~BCLO»JCLD»PCVIURAO»HOR A — L7- 

ESTABLISH  INITIAL  CONDITIONS  A 18 

2 A VERTICAL  SEPARATION  OF  THE  TWO  OBJECTS  A .19 

01  A INITIAL  SLANT  RANGE  DISTANCE  A 20 

EL  A ELEVATION  ANGLE  FROM  THE  SLOWER.  XB<  TO  THE  FASTER- XA<.  - -A--21 

OBJECT.  A 22 

.XBl.YOI.ZBI  1 INITIAL  P.QS.IUON..OF-B,, A„23. 

VA.VB  A SPEEDS  OF  A AND  B.  VA  MUST  BE  GREATER  THAN  VB.  A 24 


MOR  A RADIUS  OF  REGION  OVER  WHICH  CLOUDS  ARE  EMPLACED.  A 28 

°CVP  A PERCENT  OF  COVERAGE. .DESJRED....... „A_2I 

BCLD.TCLD  A BASE  AND  TOP  OF  CLOUD  ELEMENTS  DIVIDED  BY  100.  A 28 

4*5000. /.3.  28  . „ . ..  .... A 2.9. 

DIA4000.  A 30 

YI2/-U.L A .31. 

ELAASINXY<  A 32 

ELAEL*180./3«  1416  . . A.,  .33. 

XBIAYBIAO.  A 34 

.ZBI.A150.  _ . ...A.  .35. 

VAA750.  A 36 

JLBJ2SD. A 37 

HORA5000.  A 38 

..F.LAGAO  . . _ -JL„3.9. 

I CONTINUE  A 40 

hlAGAFLAGSl  A„„41. 

KKKAO  A 42 

ECYRt.5. A 43 

RADA4.  A 44 

.2  CONTINUE.  . . - .A.  .43. 

UCL0A4.  A 46 

-T.CLDA.9.  . - . . .K  . JtZ 

3 CONTINUE  A 48 

I.HCLQATCLCLzQCLQ.— A. 49. 

DO  6 MAI ,2  A 50 

.CHANGE  THICKNESS  .OF  CLOUD  WITHOUT.  CHANGING  .BASE  OS..P.OSIT.tUN..OF  .A.  5.1. 

ELEMENT.  ***  A 52 

rCL0A.IClDStM-K*THCLD/2.  A.  S3. 

IF  XM.GT.K  GO  TO  4 A 54 

1F~-XKKK.«GI..0.<_GU.  .10  ,4. A 55. 

CALL  SET  A 56 

USE  A NON-OVERLAPPING  CLOUD  FIELD  ONLY.  A 5.7 

DO  Not  CALL  CLOUOO  WITH  THIS  PROGRAM.  A 58 

CALL  CLOUD NO  A 59 
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4 CONTINUE  A 60 

CALL.CAICtil k 61- 

CALL  GILO  A 62 

.0Q^J«l,12. ..  A ,63.- 

A/6»30.0*tJ-l<  A 64 

— . -CALL  CLR  A 65.  . 

AZAf A28C25.0  A 66 

_C QEC1M  on  a ntRFf.T  i nt FRr.FPT  mtiNTfi<  na  Tin  rjusiwc  XHUMTr.r A 6J_ 

C SCENARIO.  A 68 

CALL.  HUNT  0 _ _ A_69_  ‘ 

KKIKK-1  A 70 

CALL.  CLKX. . A™71_ 

CALL  FREE  A 72 

AZAlAZft=2£wQ A 7 V 

CALL  HUNTO  A 74 

KMKK-1  ... A_  .75-- 

CALL  CLRX  A 76 

— . . . CALL  FREE  __ __ .A_7.7__ 

CALL  GIU  A 78 

— 5-..CONIJNUE.  ...  - — A — 7-9— 

CALL  G1L2  A 80 

CQVERAJJtXRAOtl.  0(U<A.t2^/XH0RlAZ^< A 81- 

PRINT  10,  BCLD.TCLO, RAO, COVER  A 82 

6 CONTINUE  ....  ....  A 83— 

KKMKKKC1  A 84 

-C CHANGE.  BAS  ES.  OF.  .CLUUD..W  l.TMOUT.  XHANG1NG..I HE_P.QSmQN_QF_T HE A 85_ 

C ELEMENTS.  A 86 

BCLOABCLDt.1.5 A 82— 

IF  XKKK.GT .2<  GO  TO  7 A 88  X 

....  — GO  -JO.  3 - — A— 89- 

C CHANGE  THE  RADIUS  tANO  THEREFORE  COVERAGE<  WITHOUT  CHANGING  THE  A 90 

X.  POSITION.  OF  THE-  ELEMENTS.  .THIS.  WORKS.  QNL.Y-.M1.TH.XLQUQNQ* A 91— 

7 RAD#3.  A 92 

LF— XKKK.»GI»3<-GU-IU-8 A — 50- 

GO  TO  2 A 94 

...  0.  CONTINUE.  . - - - _•  A-  .95- 

GO  TO  31, 9<,  FLAG  A 96 

9..C0NT  INUE  . — A.  .97.  - 

C A 98 

1.1, 11H  .RADIUS  » ,F5.1,13H  .COVERAGE  * ,F6.3<  A 100 

- - end  . . .. ,A..101=. 
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Table  D-3.  Listing  of  L0S1AS. 


v 


PROGRAM  LOSIAZX INPUT , OUTPUT < A 1 

C IHIS  PRUGKAK..WAS-UES  IGNED  TQ..EXAM1NE— THE— T.1ME..HLSTURY..OF— THE  .CLOUD  A 2_ 

C OBSTRUCTED  OR  CLOUD  FREE  LIME  OF  SIGHT  BETWEEN  TWO  OBJECTS  MOVING  A 3 

_£ tM_THE-PR£-SLNCF_QE  l /.LUmi  F-1F1  O— 1HE-USER.  MUST.  SPFr.lFV  THF A 4 

C INITIAL  CONDITIONS  AND  DECIDE  ON  ONE  OF  TWO  CHASE  MODELS  AND  ONE  A 5 

C OF  TWO  CLOUD  MOOELS  TO  BE  USED.  IN  THE.  ANAL  YS  IS— -AERIAL.. COVERAGE..QF  — A...  t. 


C THE  CLOUD  IS  DETERMINED  BY  t HE  FRACTION  OF  THE  UNDERLYING  SURFACE  A 7 

C.  COVERED  BY  THE  . PROJECTION  OF  THE  .CLOUD  . BASES-  - .. — 

COMMON  /SLK1/  XCLDtlOOOC, YCLDXIOOOC, JJ  A 9 

COMMON. /BLK2/  PXAXIOOC ,PYAX100<,PZ AX  100<»PXBtlOO<tl>YatlOOC»R2:aXLQO A— 10 

i<,TABX100<,0AB2100<,KK  A 11 

COMMON  /8LK3/  A2 A. AZS.EL .XBI.YB1 »Z01,0l»VA,VB  A.  12. 

COMMON  /BLK4/  KLNX200<,KL,L1 MAX4<.LIMBX4<  A 13 

-COMMON  /BLK5/ -CL  OB  . -A— .14. 

COMMON  /BLK6/  J A 15 

CQMMQM-/8LK.2/.-BCL.0.JCLU+PCVR— RAD*HOA A 16— 

C ESTABLISH  INITIAL  CONDITIONS  A 17 

C 2.  t VERTICAL  SEPARATION  OF  .THE. . I BO-OBJECTS A.  -18- 

C 01  # INITIAL  SLANT  RANGE  DISTANCE  A 19 

C EL  • ELEVAT.IQN  ANGLE. .FROM.  .THE.  SLOWER.  XB<.,IO..THE..F  ASTER. -XA<.  A .20.. 

C OBJECT.  A 21 

-C XB1-YBI.2BI— A-1N1J.1AL-PQ.S1T1QN-0F— B. A 22- 

C VA.VB  # SPEEDS  OF  A AND  B.  VA  MUST  BE  GREATER  THAN  VB.  A 23 

C-  RAD  # RADIUS  OF  CLOUD-ELEMENTSAIOO. ---  - A-  .24'— 

C HOR  # RADIUS  OF  REGION  OVER  WHICH  CLOUDS  ARE  EMPLACED.  A 25 

C PCVR  • PERCENT.  OF  COVERAGE  DESIRED.  _.A — Zb- 

C BCLD.TCLD  # BASE  ANO  TOP  OF  CLOUD  ELEMENTS/ 100.  A 27 

2# SOOO-/3. 28 A — 28- 

Dl #4000.  A 29 

SR#D  I A-  -30 

Y#Z/SR  A 31 

EL#ASINXY<  -A~,32~- 

ELIELM80./3.1416  A 33 

XUI.#.m#CL. .. A 3.4— 

281 #150.  A 35 

YAI25Q. A..  .36 

VA#1000.  A 37 

VB#.25.0.„  A 38— 

HOR#5000.  A 39 

FLAGiO-  — — A 40_ 

1 CONTINUE  A 41 

LLAliJtELAGU A 42  - 

KKKiO  A 43 

P-C.YM.,.7 .A — 44 — 

RAD#3.  A 45 

HaOJA.. . . A — .46— 

BCLO#2.5  A 47 

mm a 4a_ 

TCL0#5.  A 49 

2L.CQN.UNUF . . „ . . _A„50- 

THCLD#TCLD-BCLD  A 51 

DQ..5..M#.i.«  2.  _ A_.  .52— 

C CHANGE  THICKNESS  OF  CLOUD  WITHOUT  CHANGING  BASE  OR  POSITION  OF  A 53 

_C ELEMENTS  ..A** A 5A_ 

TCLO#TCLOGXM-l<*THCLD/2.  A 55 

IE  -XM.GU  1.<.,GU-T.U.,3 - A.  —56.  . 

IF  XKKK.GT.O<  GO  TO  3 A 57 

CALL  SEJ . ..  A-.5B- 

C OECIOE  ON  USING  AN  OVERLAPPING  CLOUD  FIELD  XCLOUDO<  OR  A A 59 

C .NON-OVERLAPPING  F1F.L  Q_tCLQULNQ£.— CLQUl)Q_iS RECQMENOED—EUR A—  5.0— 

C COVERAGES  OF  .50  OR  GREATER.  A 61 

.CALL  CLQUDQ.  - A.  62 
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3 CONTINUE  A 63 

CALL-CA.T.CUl - A- -64 

CALL  GILO  A 65 


c 

r. 

AZB#30.0*SJ-1< 

...CALLXLR,  . 

AZA#A7OC25.0 
_ .AZ.AIAIB.C155.. 

0EC10E  ON  A OIRECT  INTERCEPT 
scenario. 

SHUNT 0<  OR  TAIL  CHASING 

A 

..  A. 

A 

. ..A. 

SHUNT C<  A 

A 

pp 

67  ! 

68  ; 

69 

70 

71 
-72 

CALL  HUN  TO 

A 

73 



_ .JKKKJKK-.L  _ 

_ 

. . . _ . A 

lit 

CALL  CLRX 

A 

75 

..  .CALI FREE  ... 

_ 

...  A 

76 

AZA#AZ0-25.O 

A 

77 

A7A#A7B-lS<i. 

A.. 

7JB 

CALL  HUN TO 

A 

79 

KKUKKr.l  . - - 

A 

-80. 

CALL  CLRX 

A 

81 

.CALL  FREE 

. ..  A 

,82  1 

CALL  GIL1 

A 

b3  i 

A CONTINUE 

A 

84 

CALL  GIL2 

A 

85  * 

COVERHPXV.R 

A 

.86 

PRINT  8,  BCLD.TCLD, RAO, COVER 

A 

87 

5._X0NtINU£ 



A, 

.88  . 

KKKidKKKCI 

A 

89 

-C... 

..CHANGE.  BASES  OF  CLOUD  .W.l  THOU T CHANG I.NG...T.HE  P.QSI.T.IUN 

OF  -THE- -A 

90 

c 

ELEMENTS. 

A 

91 

nri  mmc.i  nr.i  .s 

A 

_92 — 

IF  SKKK.GI  .2<  GO  TO  6 

A 

93 

GO  10.2 

— 

..A 

9.4 

6 CONTINUE 

A 

95 

GO  TO  11,7C,  FLAG 

. 

. ..  .A 

96 

7 CONTINUE 

A 

97 

c 

. A- 

- 9H  — 

8 FORMAT  ?// , 10X, 29HAB0VE  INFO 

BASED  ON  - BASE  # ,F5.1 

, 8H, TOPS  H ,F5  A 

99 

1.1, UH  .RADIUS  H , F5-J  ,13H  , 

COVERAGE  A ,F6.3< 

_ A. 

100 

END 

A 

101- 
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_ SUBROUTINE  SET  . ...  A i 

C THIS  SUBROUTINE  INITIALIZES  REQUIRED  OPERATING  VALUES.  ***  A 2 

CALL  CLDNll A I- 

CALL  CLOOVK  A A 

CALL  HUNTS  . _ _ A 5 

CALL  FREES  A 6 

CALL  MAPS  . _ A Z 

RETURN^  A 8 

END A 9Z 


SUBROUTINE  CLONO 

B 

i 

c 

THIS  SUBROUTINE  GENERATES  THE  CLOUD  FIELD  TO  BE  EXAMINE [).♦♦***•* 

B 

2 

c 

♦CLONO  PRODUCES  A NON-OVERLAPPING  CLOUD  FIELO  AND  CAN  BE  USFD 

B 

3 

c 

FOR  COVERAGE  UP  TO  ABOUT  45  PERCENT.  SUBROUTINE  CLOOVR  IS 

B 

4 

c 

RECOMMENDED  FOR  COVERAGES  GREATER  THAN  45  PERCENT. 

B 

5 

- . COMMON  /BLK1  / XCL.Utt.OOOC. ,.YCLD'ilOOQ<.  J J 

B 

6 

COMMON  /BLK7/  8CLD.  t'CLO,  PCVR  ,R  AD.HOR 

B 

7 

c 

INITIALIZE  *«**♦ 

B 

8 

SML#0.0001 

B 

9 

JJ#0 

B 

IP 

JPSS  #2000 

B 

1 1 

PY#3.l4l6iCDR#PY/I80.0tCRD#l  .O/CDR 

B 

12 

XKA#2.*UinR-R/U)«100.<*FRA#RAO*100.iORA#2.*FRA 

0 

13 

CVR#PCVR$CL0S#0.0iBMP#PY*FRA**2 

B 

14 

TAR#PY*H0R»*2 

a 

15 

CP#SECONDf.CP< 

B 

1 6 

CALL  RANSCT  XCP< 

n 

17 

c^ 

**i.f  *><*»*:**«♦  END  OF  INITIALIZE  *»*»**♦ 

a 

.10 

RETURN 

B 

19 

ENTRY  CI.OUONO 

0 

2 f 

PRINT  7 

B 

21 

_1 

CONTINUE 

R 

.22 

c 

***  I EST  FOR  PERCENTA  LOUD  COVER  REACHED  ***** 

a 

23 

PCNI#CLDS/TAR 

0 

24 

IF  HPCNT.GT.CVIK  GO  TO  5 

B 

25 

2 

CONI INUE 

B 

26 

C 

TEST  FOR  OVER  2000  PASSES  ♦** 

B 

27 

JPSS#JPSS-I 

B.  28 

IF  TIJPSS.LT.  1<  GO  TO  5 

B 

29 

|<PSS#2000-JPSS 

B 

10 

c 

****  OBTAIN  X AND  Y COORDINATES  FROM  RANDOM  NUM  GEN  *** 

B 

31 

RX#RANFHPY<$RY#RANF?.PY< 

B 

32 

XP#RX*'XRA-tHOR-R  AD*100.<  tYP#  RY*XRA-HHOR-RAD*  100. < 

n 

33 

RP#0.0$ARG#XP**26YP**2 

14.- 

IF  HARG.GT  ,SML<  RP#SQRTHARG< 



B 

35 

C 

♦♦TEST  FQR  X AND  Y WITHIN  ACCEPTABLE  LIMITS  *** 

B 

36 

IF  5RP.GT . tHOR-RAD* 100 .<<  GO  TO  2 

B 

37 

C 

TEST  FOR  OVERLAP  ON  OTHER  CLOUDS .♦**** 

B 

30 

J#1 

B 

39 

-3 

cunuhul  . - . 

— . -11 — itjQ 

IF  ?J.GT.JJ<  GO  TO  4 

B 

41 

XJ#XCLOCJ<$YJ#YCLDHJ< 

B 

42 

XT#XP"XJ*YT#YP-YJ$RT#0.0 

ft 

43 

ARG#XT**26YT*V2 

Q 

44 

IF  HARG.GT.SML<  RT#SQRTHARG< 

ft 

4 5 

IF  HRT.LT.ORA<  GO  TO  2 

B 

4 0 

J#JC1 

B 

4 / 

....GO  TO  3 

a 

4.a„ 
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4 CONTINUE 

C . ♦**  .X..ANQ  Y .HAV.E  .BEEN.  ACCEPTED.  ***• _ 

C ♦ STORE  RESULTS  ***UPOATE  CLGUU  COVER  ***PRINT 

.JJUJJ.tl.  _ _ 

XCLOtJJ<*XPiYCLO*JJ<i»YP 

XL0S.IC  LOJS.S.aH£ 

C PRINT  8, JJ,XP,RPf CLDS,PCNT,KPSS 
GO  10  1 

5 CONTINUE 

P.CVMPCNT  

SCLfO.Ol 

HQ.  6_  JA.U.JJ. 

C * STORE  SCALEO  VALUES  OF  THE  CLOUO  COORDINATES  **'■> 

XCL0*J<#^CLDSJ<*SCL$.YCLD?J<1(YCLDSJ.<*SCL  .. 

6 CONT INUT 

PRINT  8 ? JJ.iXP.j  Y.P.t  RP«  CLDS » P.CM.T i KPSS 
RETURN 


8 49 



8 51 

a sz„ 

8 53 

B 54 

B 55 



8 57 

. _ .a  _ xa._ 

X)  59 

h 6lL_ 

8 61 

-62 

8 63 

B .64. 

B 65 

E 66 


C B 67 

C _ - .8.  -68 

7 FORMA!  Y1H1/10X,2HJJ(11X,1HX,11X,1HY,1\X,1HR,3X,4HCLDS,8X,4HPCNT,6  B 69 

1Xi4HKPSS<  _ ...  . ..  B .7.0 

8 FORMAT  «2X,I 10.4F12.  I.F12.3, 1 10<  8 ?1 


SUBROUTINE  CLPOVR 

C 

1 

c 

THIS  SUBROUTINE  PRODUCES  A CLOUD  FIELD  WHICH  ALLOWS  FOR 

SOME 

C. 

2 

c 

OVERLAPPING  OF  INDIVIDUAL  ELEMENTS.  ANY  NEWLY  INTRODUCED 

ELEMENT 

c 

3 

_c... 

I S_  ALL  OW  ED.  T O _0V  ERLAP..NQ  OTHER  .OR.  ONE  QT  HER.CLQL)  DELEMENT 

. IF.  THF 

r. 

4 

c 

NEWLY  INTRODUCED  ELEMENT  WOULD  OVERLAP  TWO  OR  MORE  EXISTING  ELEMEN 

c 

5 

c 

IT  IS  REJECTED.  THE  RESULT  IS  THAT  STRINGS  OF  CLOUDS  ARE 

POSSIBLE 

C 

6 

C 

BUT  GROUPS  OF  MORE  THAN  TWO  CLOUOS  ARE  NOT  POSSIBLE.  *** 

c 

7 

COMMON  /BLK1  / XCL0*1 000<  . YCLD'UQOOC,  J J 

c 

8 

COMMON  /BLK?/  OCLO, TCLO, PCVR ,RAO,HOR 

c 

9 

♦ **<.**IN!T  ITIAI  1 7P  +♦*«¥*  

c 

l n 

SML  H 0. 0001 

C 

II 

JJ#0 

C 

12 

JPSS/WPSS#  1000 

c 

13 

PY#3.1416$C0K#PY/180.0$CRU*1.0/£UR 

c 

14 

XRA/H2- ♦-"HOR-RAD*  100. < 

C 

15 

- r-„ 

ERAtfRAQMQQ.. . _ . 

_c 

_L6_ 

URA  # 2. * FRA 

c 

17 

CVRtfPCVR  . - - 

c 

18 

BMPUPY+FRA*,' . 

c 

19 

rARm*H0R**2 

c 

20 

CLDS0O.O 

c 

21 

- - .CfiHSEC0ND7CP.<  „ ... 

■ — - 

_c_ 

-12- 

CALL  RANSET  *CP< 

c 

23 

c 

«.+*«.*#.  + •***»*  END  OF  INITIALIZE  +•****»* 

c 

24 

RETURN 

c 

25 

ENTRY  CLOU DO 

c 

26 

PRINT  9 

c 

27 

1 

CONT  INUF 

c 

23 

C 

Hi  TEST  TOP  PRRCENT AGE  CLOUD  COVER  REACHED.  1((« 

c 

?9 



PCNT  #CLl)S/TAR  

___  - — i — 

JL. 

-IS- 

ih  *pcnt.gt.cvr<  go  to  6 c ji 
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CONTINUE  C . 32 

**  TEST  FOR  OVER  1000  PASSES  ****  C 33 

JPSSflJPSS-1  _ C 3A_ 

IF  2JPSS.LT. 1<  GO  TO  6 C 35 

KPSSIMPSS-JPSS  C 36 

****  OBTAIN  X AND  Y COORDINATES  FROM  RANDOM  NUM  GEN  ***  C 37 

RXflRANF2PV<$RYtfRANFSPY< 1 C 38 

~XP#RX*XR  A-2H0R-R  AD*  100..<  C 39 

YP#RYJXR_A-2HOR-RAD*100.< _ C AO 

RP#0~.0SARG#XP**2£YP**2  C 41 

IF  2AKG.GT  »SML<  RP ASQRT2 ARG< C 42 

**TEST  FOR  X AND  Y WITHIN  ACCEPTABLE  LIMITS  ***  C 43 

IF  2RP.GT.  2H0R-R  A0*JL00 . <<  GO  TO  2 C 44 

*TEST  FOR  OVERLAP  ON  OTHER  CLOUDS.  ***  C 45 

J#l_  C 46 

KBMP# l C 47 

BMPXflBM  P C 48 

CONT INUE  C 49 

J.F_  2J.GT.JJ<  GO  JO  5 C 50 

XJtfXCLDi j<JYJ#YCLD2 J<  C 51 

XTtf XP-X J $YTtfVP-YJ$RT#0.0 __ C 52 

ARG#XT**2£YT**2  C 53 

IF  2ARG.GT.SML<  RT SQRT2ARG< C 54 

IF  2RT.LT. 0RA<  GO  TO  8 C 55 

CONTINUE  _ C _56 

JtfJCl  C 57 

-GOJLO  3__  _ C ,58 

CONTINUE  C 59 

»*«  X AND  Y HAVE  BEEN  ACCEPTED  »*♦* C 60 

C * STORE  RESULTS  ***UPDATE  CLOUD  COVER  ***PRINT  **  C 61 

_ __JJ#JJJ1  _ C__62_ 

XCL02JJ<#XP$YCLD2JJ<#YP  C 63 

CLDS  tfCLOSSBMPX  _ _ C 64 

GO  TO  1 C 65 

6 CONTINUE C 66 

PCVR0PCNT  C 67 

SCL  #.0.01  . C — 6.8_ 

DO  7 J#1,JJ  C 69 

_JLCLD2J<S  XC  LD2J<*SCL$YCL02J<fl  YCL02J<*SCL C 70 

7 CONTINUE  C 71 

PRINT  10,  JJ.XP.YP.RP.CLDS.PCNT.KPSS C 12- 

RETURN  d 72L. 
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CONI  1NUE 

U'Sl  I- OR  UVTRLAP  ON  MORI;  1HAN  (INC  l LEMON  I • REJECI  IF  I HI5  NtW 
ELEMGN1  OVERLAPS  MORE  THAN  ONE  OlllLR  tLEMENl  , OR  IF  II  OVERLAPS 
RY  A UIS1ANCE  ORE  AI  CR  1 MAN  1 HE  RADIUS  OF  T HE  OVERLAPPED  tLEMENl.  ♦* 

_ilxrl.li,,fm<.  go  lo  2 . 

IP  7KRNP.L I . l<  00  HI  2 

‘DETERMINE  I lit  ADDITIONAL  COVERAGE  INDUCED  0 Y HIE  P01LNUALLY 
ACCEP1 A RLE  NEW  ELEMENT « NOl F THAI  IMF  ADDITIONAL  COVERAGE  WILL  NO 
RE  USED  II  A MOVE  IS  MADE  10  ADDRESS  2 BEFORE  A MOVE  IS  MADE  10 
ADDRESS  5.  *** 

01HRT 

02HDRA-IU 

_.UJ#£RA-D2/.2,0  _ . . - 

DA#  0.  Ot  A RGKERAt  * 2-03“  ?. 

II*  XARG.Gl.SML<  DA#SOKUARG< 

ARG#0A/0RA 

ARC  If  DA  / OR  A TANG  # 2 . 0 * A S 1 NX  ARCS 
RAT  # ANG/T2  »0*PY < 

-Al#  D3AUA/.2..  0 » A2KRAI*.  RMIU  A J#  2..  Q>.XA  2- A L<.  ...... 

RHPXfUUMPX-A) 

RUMP  Ilf  RUMP- 1 
GO  10  A 


FORM  A l UKl/  10X.2IUJ,  ttX,  1UX  , 1 1 X , IHY , It  X , 1MK ,8X. AHCLOS, OX, AHPCNI , 6 
l X,  AHKPSS< 

FORMAT  T.2X.I  10, AF 12. I » E 1 2 .3,  110< 

END 


C 

C 

C 

C 

X_ 

C 

C 

C 

C 

C 

C 

c 

£ 

C 

C 

C 

C 

C 

C. 

C 

C 

c 

C 

C 

c 

c 

c 

c 

c 


7A 
lr> 
lb 
77 
_.7_ U. 

7 V 
DO 
in 
02 
in 

RA 
OS 
. «6„ 

8 I 
OR 

09 

go 

91 
1)1 
93 
9 A 

•r> 

96 
9 1 
.9 11. 
99 
100 
10 1 
10?- 


.1  111:  VEHICLES  IS 


SURROUTINE  IIUNl  F 

C * THIS  SUBROUTINE  DETERMINES  1 HE  IH3R LIONS  UI- 
C A FUNCTION  OF  SIMULATED  TIME. 

COMMON  /RLK2/  PXAX 1 00< ,P YAX100< . PiAil 00< ,P XRXIOO <, P YRX 100< , P 7 11X100 
1 <,TARX100<  *DARX100<»KK 

XOMMON  VJILIO/..  A2A. AIR , EL. XR L.XU l..IRl.,DL,VA ,.V.R 

COMMON  /U l.KA/  KLNIL’OOC.KL.LI  MA*4<,LIMRXA< 

COMMON  / RLK'i/  CL  Dll 
COMMON  /BLK6/  J 

COMMON  / IK  7/  I1CLD.  I OLD,  PCVR  .KAO.HOR 
C INI  V 1 AL I Li . ***** 

. „_$MLf  0.0  0001  . 

PY# 3 . 1 A 1 61  COR#  PY  /l  RO.OiCRDffl  .O/COR 
C END  OF  INI  1 1 AL1 2E. **♦+* 

RETURN 
ENTRY  MUNI 0 


0 

D 

IT 

0 

IT 

D 

0 

0 

D 

D 

D 

D 

0 

0 

0 

D 


1 

2 

1 

A 

r> 

6 

r 

R 
9 
10 
1 1 
12 
1 3 
1A 
l’> 
16 
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C DETERMINE  POSITIONS  FOR  DIRECT  INTERCEPT  MODEL  USING  T IMF  STEPS  OF  0 17 

C_  .0.2— SECONUS NUTE T.HA.I 1 N_1  IUi_.MUDET__IIiE_EL£UAJJUN_ANGLE_REMAIJJS D__LH„ 


C CONSTANT.  **♦  0 19 

Pill  #EL  0 20 

C PRINT  13, AZA,AZU,GL,01,ZBI  0 21 

A*C0RtAZlHVXl3llCUStA<»Vl\JVYll#SlNtA<*VU  0 22 

A#COR*AZAiEi*COR*EL  0 23 

_ JCAI.AUI.tCQSXLCt.COSJAC  — 0-.2-A- 

YA1  *DHCOS'tl.<»SlNtA<  I)  29 

ZAI#0I<,S1NXQ<CZUI  D 26 


K XA#XAISYA#YAliZA#ZAl  I)  27 


• XBfcXUIiYBfcYUItZU/KZUIilHIUI  0 26 

CLDDrfBCLO  0 29 

.KPSS#0 — D 10- 

< T MiKO.O  0 Jl 

JTMUOiJOLI #20  0 32 

DT#J0LJ/lO0.0  0 33 

KK*0  D 3A 

C PRIN1  It  D 39 

KF.LAGIO  0 36 

KMKL/90C.1  0 17 

_C ULLERtlLNE  me.  hi  N A I X.tX-EU  S.l.U.UN-  Oh  JLtUL  .VE1U.CLLS..S  Si. 0.  _XU__ 

A#CDR*AZtl  0 39 

DUM 1 #C0StA<*  XXA&YA''  1 AN'S  A <<  D AO 

I)UM2WA**2/VB**2<-1.  0 A 1 

UB«.*-OUMlf.SQRIXOUMl**2f.UUM2'‘Dt*<2«/OUM2  U A2 

xi  #ou*cns7:A<  o a j 

XfcJLQU*$JMtA< CL _AS.. 

OA#  SORT  t T.XF-XA<**2CtYF-YA<+*  2CXZI3-ZA<*'12<  0 A9 

PUXr-XA</0A»0«HYP-YA</0Al.R#tZU-ZA</UA  D A6 

1 CONI  INUF  0 A 7 

C DETERMINE  POSH  IONS  OF  A AND  II,  FOR  THE  FINAL  POSH  ION  JUST  CALCUL  0 AU 

C A TED,  EVERY  0.2  STCONUS. *****  0 A9 

KMfKKG  1 _ „ 0 90 

PXA  TKK<rf  XA>PYAXKK<  HY  AtPZ ATKK  <#Z A 0 91 

PX0'tKK<.YXl'VPYI3'£KK<#Yl3$PZD1KK<#Zl3  1)  92 

T A»-«KK<#TM  HlAl3«KK<#t)  n '.1 

C PRINT  U.KI’SS.TM.XA,  YA,ZA,X0,Y6,Z6,D,PMI  U 9 A 

IF  *KI LAG.LO. 1<  GO  TO  2 0 9*, 

-C 1F..A.  .HAS.  RLACIltU.  EASE.  .OI‘..CLQUUS,..tLAG-.TIl.LS.JLOSlIJiQN  .1).  9.6. 

IF  tZA.LT .CL  1)13+  100.<  GO  TO  3 0 9 7 

2 CONTINUE  I)  9}| 

KPSS#KPSSf.  1 n 99 

IF  XKPSS.GT . 100<  GO  TO  A U 60 

C IF  SLANI  RANGE  IS  LFSS  I HAN  901  ONUS,  SlOP.*****  0 M 

.C  SAV.L  LIU  S 1’US.U  ION .****<  . . . 0.  62 

IF  'll). LI  .901  .0<  GO  10  A 0 6' 

JTMdJTMCJDU  HMUJTMTTMinM/lOO.O  U OA 

IF  t IM.G T . 9. 0<  GO  TO  A 0 69 

VXA#VA*P»VYA#VAMjmAm*R  L'  66 

XA#XAf.VXA*DT$YAIfYA£VYA*I)It/A*/Af.  V7A*D1  0 67 

Xm04YXl)*Dl»YR.m6VYU*0mi»#ZU  0 6.6 

l>*0.0  0 69 

ARG#XX0-XA<**2C  iYI'-YA<  **  2LTZl'-ZA<**2  D 70 

IF  t AKG.GI  .SMI. < 0# SORT tARG<  0 71 

GO  10  1 0 72 

3 L lMlltKM<#KK-l  0 7 1 

KFLAGUl  II  7 A 

GO  10  2 D 79 
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A IP  XL IMAXKMC.LE. l<  LIMAXKM<IKK 
IF  XLIMtlXKM<.LE.  1<  L 1MUXKM<#KK 

scmo.oi 

C SCALE  THE  POSH  ION  VALUES.  ♦** 

DO.  ii_Kli..KX __  ___ __ 

pxaxk<ipxaxk<*scl»pya*k<ipya*k<*scl 

PZAXK<*PZAXK<*SCL»PXBXK<#PXUXK.<*SCL. 

PYBXK<l»PYHS!K<*SCLiPZBXK<#PZUtK<*SCL 

oabxk<#oabtk<*scl 

5 CONTINUE 

RETURN 

ENTRY  HUN  I 

C IMIS  ENTRY  DETERMINES  THE  POSITIONS  OP  1 HE  I HO  VEHICLES  WHEN  A 
C ALWAYS  HEADS  TOWARD  THE  CURRENT  POSITION  OF  B. 

C PRINT  13  • AZA  , AZB,EL,01,ZUl 
CLUB0HCLD 

AHC11R*  AZ  111  VXBUCUSt  A<+VUiV  YU  < SINXA9YVU 
A»COR*AZA*E»CI)R*EL 

XAIJ DISCUS XE.<1CUSXA.<.  __ 

YAI  #1)1  *COSXE<*Sl  NXA< 

ZAIiOltS INXE  < 

XAdXAI tYA#YAliZA#ZAl 
XU*  XBI  iYBKYB  I IZtiKZBI  iUKU  I 
KPSS  #0 

PHI# EL - 

TM#0.0 

J1M#0JJ0L1#20 
D I # JOLT / 100. 0 
KK#0 

C PRINT  11 

KELAGHLELAG#  Q.  . ....  . . . 

C ESTABLISH  THE  BOUNDS  WITHIN  WHICH  STATISTICS  WILL  BE  COLLECTED. 
UPPERKELCS. 

LOWERHEL-9. 

KMIKL/5QC l 

6 CONTINUE 

_C. DETERMINE  .1  HE  POSITIONS  .OP  A .ANU.  B.  _ 

KK#KKf.l 

PXAXKK<#  XAJPYAXKKCIIY AiPZAXKKCAZA 
PXU  IKK<#XIHP  YI1%KK<  #YI1  JPZ  BTKK<  #ZB 
TAB  v.KK<«TMi0AUXKK<#0 
IP  XKPLAG.L0. l<  GO  TO  7 


l)  76 
0 77 

o / a 
I)  79 

U R0„ 

0 HI 
D (12 
D 8 3 
0 BA 
0 89 

U 86.. 

D 61 
0 Bit 
D 89 
0 90 

0 91 

U 92 
0 93 

U__9.4._ 

l)  99 
1)  96 

0 97 

D 98 
I)  99 

tL_100._ 

0 101 
U 102 
I)  103 

n ioa 

0 109 

U_LQ6. 

*'  I)  107 
1)  108 
l)  109 
D 110 
Dill 
.11 -LL2 
D 113 
U 1 1 A 
D 119 
D 116 
OUT 


-C.  IE  A HAS  REACHED  THE  -BASE  .01  JHC  .CLOUDS,.  FLAG  THIS  P-OSI  UUN.  0.118 

IP  IZA. LT .CLUB*  l 00 .<  Gil  TO  8 I)  119 

7 CONTINUE  0 120 

C PRINT  12,KPSS,!M,XA,YA, ZA,XB,YB,ZB,D,PHI  U 121 

KPSS#KPSSC1  0 122 

C IP  THE  SLANT  RANGE  DISTANCE  IS  LESS  THAN  901  UNI  IS, OR  IP  TOT  D 1 2 1 

C SIMULATED  1 1 ME  HAS  EXCbtOLll  9 SECONDS,  SHIP.  Ml  0 129 

IP  XKPSS.G1 .90S  GO  10  9 D 129 

It  XIMI.90l.tK  GO  TO  9 D l .» <. 

jlM#JIMr.JDLTilM#JTMUM#lN/lOO.O  D 121 

II  JIM. GT. 9. 09  GO  TO  9 0 128 

P#JXB~XA</DEU*tYB-YA</IHR*17  B-  ZA</P  D 12‘» 

VXA#VA*PiVYA#VA*01V2Af!VA*K  U 1 1U 


XAl*XAf.VXA*01  1-YAH  YAGVYA*H HZA «7  ACV/ A<*D1  0 Ml 

XBK  XIIT.V  XB*D  UYB  YUCV  YU  *0T  TZiUtZB  D 1(2 
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C 

C 


DETERMINE  IHE  CURRENT  ELEVATION  ANGLE.  *** 

0*0.0 

ARG*XXB-XA<**2£XY8-YA<**2£XZB-ZA<**2 
_IE — XARG.GI..SMLC-UXSORlJtARG< 


SCALE  THE  POSITION  VALUES. 

SCLJIO.O.l 
OU  10  K*1»KK 

PXATK<#PXAXK<«SCLiPY6XK<#PYAXK<*SCL 

PZA*K<*PZAXK<'*SCLiPXt)'tk<»PXntK<*SCL 

CL  m BJUiltRZ  i)J.lL<±SC.L 

OABXK<*l)ABXK<*SCL 

CONTINUE 

RETURN 


-tNQ 


D 133 
0 134 
D 135 
-0-136- 


PHI IATAN2XK.  SDK  TXP**  2CQ**  2«  D 137 

PHI  *-PHI  *CRD  -0  -138 

IF  XLFLAG.EO.K  GO  TO  6 D 139 

C IF  THE  ELEVATION  ANGLE  HAS  EXCEEDED  THE  ANGULAR  LIMITS  ESTABLISHED  0 140 

C FLAG  IMIS  POSUtON.  ♦**  0 141 

1E-XPH1-.LT-.  UPPER.  A.  P HI.  GT  .LOWERS— £Q— TO— 6 0-143- 

LIMAXKM<*KK  0 143 

LFLAGA1  0 144 

GU  TO  6 D 145 

8 LIM8XKM<*KK-1  0 146 

KFLAGI1  0 147 

GO  TO  7 ...  D 148. 

9 IF  XlIMAXKM<.LG.l<  L I MAXKM<*KK  0 149 

IJE_SUMi)Jl!Ji-<.LC..l<.  .LIUIIS1SMSAKK. D_i£Q_ 


I)  151 
Q .1.52.. 
0 153  , 
.0.  15.4.  , 
0 155  1 
JL15.6.. 


I)  157 
0 15.0 
D 159 
D 160 
D 161 
-D-162=. 


SUBROUTINE  FRFEF  E 1 

C THIS  SUBROUTINE  USES  THE  POSITIONAL  INFORNAI  ION  PRODUCED  IN  .El.THLR  E 2 

C HUNT  OR  HUNIO  AND  THE  CLOUD  INFORMATION  GENERATED  BY  EITHER  CLDNO  E 3 

_C  - SB.  CUIOB.  10  .DEI ERMINE.  THL.Jiil.SJ.UJiX  -OF  .111 E._C£LQS-B ETJrt EEN_A-&E10LJI_^_! — E 4_ 

COMMON  /Bl  K l / XCLUX1000<  »YCLDX1000<» JJ  E 5 

COMMON  /BIK2/  PXAX  100<  »P  YAX100<*  PZ  AX  100<  tPXBXIL'O  <»PY.Bt  100<»  PZBX100  E 6 

l <*  T ABX100<  tOABX 1 00< , KK  F 7 

COMMON  /BLK4/  KLNX200C. KL,LI MA?4< ,L l HBX4<  E 8 

COMMON  /BLK5/  CLOU  E 9 

. . _ -Cl!MMQ£L7J)LJ$<i/  J _ - E_1.IL 

COMMON  /0LK7/  BCLD, TCLO, PCVR ,RAD,HOK  E 11 

C INITIALIZE-  **♦  E 12 

PY*3.l4I6*CQR*PY/lflO.O*CRO*l ,0/CDR  C 13 

SML# 0.0000 l C 14 

CR#R  AUtCZV  2.*CR  E 15 

SMLHA.0..01  ...  E — Lii_ 

KBLNKV 1R  c 17 

KX* IRX  E 18 

C END  OF  INITIALIZE.  ***  E 19 

RETURN  C 20 

ENTRY  FREE  F 21 

- CLDBJBCLDKLOWCLO  - . _ ~ .-JL.-22 

JFLG*KFLG*0  F 23 

CR»RAU»CZ#2.*CR  L 24 

KNT  *0  F 25 

h»0  E 26 
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CONI  I NUP  0 27 

KNUlv.NI Cl  _ . _ __ E—21L 

JCrfl  b 2 V 

******  SUUl  UK  Hi  IN  TU  AX,  AY  AND  ROT  AIL  AXES  PUR.  11Y#0  ****  t 3Q 

SLOPE AH/D  *****  c 31 

KrfKCl  E.  32 

ib  £K.GI.KK<  0.0  TO  17  E 33 

JvXft  PXAXKO  A Y ‘.M  f’ Y A XK.O  ALYL#  P.2AY.K.C-  _ E 3A_ 

BXAPXBtK  < 1 BY  APYB  'SK<P0LVL  APZB'SK<  E 15 

HAU Y-AY 3 0XBX— AX VRASOR r£li*"2f.O*'  ?<  E 36 

SGNA1.0UI  TO.l  1 ,0<SUN*-1.0  E 37 

lb  SABSSIX.LT.SMLIK  l)#SGN*SMI.II  E 3B 

Si.  PAH/D  E 39 

SI.  PZ  ft  f.GL  VI.  - ALVL  </R  E AO 

ANG*  AIAN27.il,  OOOGRSAANG*CRO  P A 1 

-LSNAXCQSiANoCir  SNAAS  INilANGC  - _ _ E _A2~ 

*♦,*♦»***♦»  Rill  ATI  UN  IL'UAIIONS  *****  E A3 

XPA  tY-AY  <*PSNAC‘JX-AX<*CSNA  ItVHtM  L AA 

YI’A  tY-AY  <*  CSNA-Y.X-AX  <♦  (SNA  E AS 

***♦  +»*'**',♦<♦*!.♦  *«»»**<i^*v\**  12  tfh 

PRINT  26  E A 7 

KNKKN1-.1.  . . ,G ,A8_ 

PRINT  2‘i  , KN , T ABY.K<  , PXA7K< , PY  A£K<  , P7  A£K<  , PXBY. K< , P YB  £K< , P Z B£K< , DAB'i  P A9 

IR<»K  E 50 

TEST  TOR  A AND  B ABOVE  UR  BELOW  CLOUD  PIELO  E 51 

IP  SALVL.GT.CL01K  GO  TO  2 X 52 

ip  •xblvl.gt.u  im<  go  io  2 E 53 

A.JVNO.,t\  .AR.t  BuTU.  BELOW  .CLOUDS P 5* 

GO  TO  11  P 55 

CONT 1NUE  E 56 

IT  tAl.VL.1  T.CLDK  GO  TO  3 E 57 

IT  SBl.VL  .LI  . CLOT  < GO  10  3 E 58 

A ANO  B ARE  BOTH  ABOVE  CLOUDS  E 59 

.GO  T.U,  IX  - C — 60. 

CONTINUE  E 61 

CONTINUE  - L 62 

CX*XCU)t:jl.<iCYlYCLO'CJC<  E 63 

1RMUCY-AYJ.1KM2ACX-AX  E 6A 

XATRMl*ESNAt'.TRM2*CSNA  P 65 

„YHTKMHCSNA.-TKM2*PSNA  . ...  .£.  .66. 

TE.sT  PUR  UO  lit  TWPEN  A AND  B E 67 

II  £X.L1.-CR<  GU  TO  6 E 6B 

IP  TX  uV.T . YRf.CR  < < GO  TO  6 E 69 

CURIO  IS  I OCA  I LD  BtIWPl'N  A AND  B E 70 

***  TISl  I OK  A INS  I OP  CYL  HASP  E 11 

RrxO.Q  ....  ,E__  .7.2 

ARl.«X*  + :*CY** r 73 

! I CARl.-Cl  .SMI  < lUASURUAKGS  b 7A 

IE  SRI. 11 .CR<  GO  TO  12  l 75 

» ?l  ST  EUR  U INSIDE  CYL  BASl  E 76 

XHX-<UYTAYUU#0.0  T 77 

ARC  AX  1 ♦ *-2CY  I ♦ *2  X 741 

IE  '{AKG.G1  .SMI  s R I ASQR  I'l  ARl.V.  P 79 

IE  SRI .LT.CKS  GO  10  l I T 80 

11  SI  I UR  l OS  lNIturiPT  OP  CLO  PRJUN  ON  UU  XY  PI  ANl  I.  Ill 

AUY  A AtvSlYv,  E 82 

IF  SAUY.l  1 .1  IU  GO  Is  7 E II  1 

CGNUNUt  X 8A 

JvAJU.l  T 8*> 

II  tJL.LI  .JJs  1.0  IDA  L 86 

i.V  10  11  I 87 
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* 

* 


6 CONTINUE 

C LOS  NOT  OBSCURED  BY  CURRENT  CLOUD 

GO -TO. -5 


7 

C 

C 

C 


CONTINUE 

LOS  INTERCEPT  OF  CLO  PRJCTN  ON  THE  XY  PLANE 
TEST  FOR  LOS  CLO  INTERCEPT  ON  THE  ZX  PLANE 
PRINT  22, JC 

STP#0.0$ARG#CR**2-Y**2 
. If  ? A K CL«.GI  • .S M L < S.LP#SQRT?ARG<  . _ , 

X1#X~STP$X2MX£STP 

Z11ALVL£SLP2»X HZ2«ALVL£SLPZ»X2 

CONTINUE 


JL£_ 


.J.l#BLVL=iSLPJj1.a2.- 
GO  TO  8 

-CONTINUE 


IF  %JFLG.GT.0<  GO  TO  15 

JF  2DAB2K-<.L.T.5.L<-G0  .TO  18 

15  CONTINUE 

- LF  ?KF.LG.GT.0<  GO  TO.  16 
IF  *DAB?K<.LT,5.1<  GO  TO 


-E— 88. 

E 89 

_£ — ao_ 


19 


GO  TO  l 


E 91 
E 92 
E 93 
-E — 94~ 
E 95 

_E 9j6_ 

E 97 
E .9.8. . 


E 99 


c 

IF  ?Z1.GT.CLD8<  GO  TO  9 
IF.  U2 .S.T .CL DB.<  GO  .TO  9 

E 101 

Z1  AND  12  ARE  BOTH  LESS  THAN  CLD  BASE 

GO  TO  5 

E 103 

F 104 

9 

CONTINUE 

E 105 

.C 

. _ . TE5  T..FOR  ,.z  1 A.ND  12  B 01 H .GP.5  A T.ER.  AH  AN.  XL.OUD.J.OP^ 

**■* 

E 106., 

IF  ?Z1.LT.CLDT<  GO  TO  10 

E 107 

. . If  .a5Z2vLT.CL01<:  G0..T.0  .10  ..... 

F 10B 

C 

Z1  ANO  12  ARE  BOTH  GREATER  THAN  CLOUD  TOP 

E 109 

GO  TO  5 

E 110 

10 

CONTINUE 

E 111 

_c. 

LQS_.LS_Q_8SlCUR£D.  * * *_.0B  J A.I  N.  JieXJUSfJ  OF_A_ANQ_JB_ 

e ii2 

KLN%KL£K<#KBLNK 

E 113 

JL. 

PRINT  24, JC 

E_JJL4 

GO  TO  14 

E 115 

LL 

CONT  INUE 

E 116 

C 

CLOUO  FREE  LINE  OF  SIGHT  FOR  CURRENT  A AND  B 

E 117 

PRINT  23j JC-  ...  . 

P lift 

GO  TO  14 

E 119 

-12 

....cqnt.inuf ...  ....  . 

E 120 

C 

PT  A FALLS  WITHIN  CLD  PRJTN  ON  THE  XY  PLANE 

E 121 

c. 

PRINT  20, JC  . 

F 177 

c 

TEST  FOR  LOS  CLD  INTERCEPT  ON  THE  ZX  PLANE 

E 123 

. .Z1MLVI 

F 174 

STP#0.0$ARG#CR**2-Y**2 

E 125 

IF— ?AR.G..  G.T.,.SML  <.  _ STP.  K SO  RT  S ARG  < 

F 1?  h 

X2WXCSTP 

E 127 

Z2#ALVL£SLPZ*X2 

F 128 

GO  TO  8 

E 129 

-13 

—.JCMUNUt.  . ..  ....  . . 

E_L3Q_ 

C 

PT  B FALLS  WITHIN  CLO  PRJTN  ON  THE  XT  PLANE 

E 131 

_Q_ 

P.RI.NL..2 l,.J.C . ..  . 

_i_132_ 

C 

TEST  FOR  LOS  CLO  INTERCEPT  ON  THE  ZX  PLANE 

E 133 

Z2#BLVL*XP#X~RSYP#Y 

F 134 

STP#0.0$ARG#CR**2-Y**2 

E 135 

L£_*ARG.G.'USML_<  S.T.P.#^QRTJARG.< 

E.-136. 

B2#STP-XP 

E 137 

E_13B. 

E 139 
F...14Q, 


E 141 
-EJL42. 


E 143 
14.4. 
E 145 
...■E....14A. 


E 147 
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.17 

CONT 1NUE 

E Htl 

KLIKLC50 

E 149 

RE1URN 

E 150 

18 

CONTINUE 

E 151 

— 

JELGtf 1 

E 152 

KENXKl£K<»KX 

E 15J 

GLL.1.Q..L5 

C 1 S 4. 

n 

CONI INUE 

t 18  5 

KF.LG#1 

E I5t> 

KLNXKl.f.K<*KX 

E 187 

GO  TO  1& 

E 158 

c 

r 1 59 

_c 

— - — - __ . 

_ „ __ 

- - - E 1 6(1. . 

ENO 

e 1 1»  i— 

SII8 ROUT  INE  MAPF 

r 

1 

c 

THIS  SU8ROUT INE  IS  USED  TO  G0LLFC1  CFLOS  STATISllCS. 

♦ *♦ 

C 

> 

«. 

-COMMON  /IUK1Z.  XCLOiinQnC.XCLnilOOQC.JJ- 

K 

. _1 

COMMON  / 0 1. K 2 / I’XAX  100<  »P  YAXl  00<*  PZA'I  100<  *PXI\t  100<*PYDTI  00<,  P/Bt  l 00 

P 

4 

.1  <*T  AUX1Q0<  ,DABX10Q<«  KK 

r 

5 

COMMON  /BLK3/  AZA,AZB,F.L 

F 

<> 

COMMON  /ULK4/  KLNX200< ,Kt  ,L I MAX4<, LI MBX4< 

r- 

7 

COMMON  / BLK6/  J 

p 

8 

COMMON  / Ml  K 7 / IlCL  llt  TCI  Ot  POUR  rR  AO  rHOU  ........ . 

_F 

. . 9. 

COMMON  /8LK8/  1,1. 1ST,K»KIIX12»4< 

E 

10 

COMMON  / 8LKK7  KK1 ,KK2,KK3,KK4 

r- 

11 

DIMENSION  KOUNPX12< 

E 

12 

K0SH«1R- 

p 

U 

c 

INITIALIZE  *** 

r- 

1 4 

.KBLNKAl R . . 

,„E. 

1 5 

KXIIRX 

F 

16 

c 

END  OF  I N 1 1 IAL1ZE-  t** 

r- 

17 

RETURN 

F 

ID 

c 

BLANK  THE  CPLOS  STORAGE  ARRAY  AND  RETURN  THE  ANGULAR 

AND 

r 

19 

c 

CLOUO  BASE  POSITION  FLAGS  TO  ZERO. 

F 

20 

ENIRY..CLR.  ...  - . 

. .«  

J£. 

. .21.- 

KL*0 

F 

22 

DO  l 1*1,200 

F 

2 3 

KLN*K*KULNK 

F 

24 

1 

CONT  INUE 

F 

25 

DO  2 K# 1 , 4 

F 

26 

-2. 

L1MAXK<*L1MB7K<*0 

- - 

-F- 

— 2-7. 

RETURN 

F 

28 

ENTRY  CLRX 

F 

29 

C 

MU  ME  l.PLUS  STORAGE  FIELD  Mill!  DASHES.  INSERT  AN 

X PRIN1  AT 

P 

30 

C 

1 HE  ENO  OP  1 HE  FIELD. 

P 

31 

DO  3 1*1, KK 

r- 

■'2 

— 

.KLNXKLUXAKUSM. 

_ 

r 

J3 

3 

CONT 1NUF 

F 

34 

MKLCKKf.1 

F 

35 

Kl.NXK<#KX 

F 

16 

RETURN 

P 

3/ 

ENTRY  GILO 

lo 

38 

C- 

ZERO  THE  CILUS  10IAL-HU-S  MILD..  *** 

F 

39 

DO  4 L*  l , 12 

P 

4 0 

4 

KIXINTXL<*0 

F 

41 

RETURN 

I 

4 3 

ENTRY  Gill 

1 

4 1 

4.' 
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C COLLECT  THE  CFLOS  STATISTICS.  F 44 

JKAIANTAIO  F 45. 

L0J  F 46 

UU-&-KI.1 . 2 E 4.7_ 

LTS  TOO  F 48 

KLMK-.1<*50  F 4.9 

DO  6 101. KK  F 50 

AI.S.LflM.Ntl-£KL<  F 51. 

IF  IKTST.NF.KHLNK<  GO  TO  5 F 52 

LLS.LIQ - ... F . 53. 

GO  TO  6 F 64 

COLLECT  TIME-RELATED  CFLOS  STATISTICS.  KNT  A 1$  USED  TO  COLLECT  F 55 

INSTANT ANFOUS  HITS.  L T S T IS  USED  FOR  TIME-RELATED  VALUES.  •*«  F 56 

KNTA0KNTAC1  F 57 

LTSTOLTSTf.  1 F 58 

I F<UST.GE.5<  S AU  „CAT.OJiL £_  .5.9. 

CONTINUE  F 60 

KOUNT  ?L<#KNTA  F 61 

RETURN  F 62 

ENTRY  G1L2  F 6.1 

GATHER  INSTANTANEOUS  CTLOS  STATISTICS.  **♦  F 64 

LSUMI.Q. FL_6A. 

KSUM0KK  F 66 

UU  7 LI l. 12  F 67 

LSUMILSUMf.KOUN  TtL<  F 68 

STATK«XLSUMU.0</tKSUM*24.<  T 69 

PRINT  8.  S1ATK  F 70 

.GATHER  J.HL....UttL-RCLA.T.LO.  .CFLI1S.  .S.LALI  S.T.l.C S .„*?*_ F„.  XL. 

CALL  GILCAI3  F 72 

RETURN  F 73 

F 74 

F 75 

F 76 

-JL.  .JLUBMAI  J/i.lOX.,  X/MPBQ.U AD1L UY  UF  A .CHLO.S.  - SLAU.C.  .CAS  E.-,.E9...1.<___  X.  7.7. 

END  F 7 8- 


SUOROUI INF  CATCH  G 1 

C THIS  SUGROUI INF  IS  DFSIGNED  TO  GATHER  THE  TIME-RELATED  CFLOS  G 2 

_.<L  .SI  A1.1S.LICS..  • **.  . ....  ft.  2 

COMMON  /RLK3/  A7A  , A7H , EL ,XR l , YU  I , 281 , 01 , VA , VB  G 4 

CUHMON  / 11LK f f llCl.O.TCLO,PCVR,RAD,HOR  G 5 

COMMON  /RLKll/  1. . L I S 1 , K, KLUI  2 , 4 < G 6 

CUMMON  /ULR2/  I'XAt  1GQ<  ,P  YAT100<  , PiAT  100<  .I’XIIXIOD  <,  PY  tit  1 0lK.P2tit  100  G 7 
l < » TAHX 100<  (DAIT'll  00< , KK  G 8 

DIMENSION  05X12..4<.  U10XL2,4<.,  I1151U,.4<,  l)2QU2.4<,  »25X12,.4<..  0.1  .ft  9 
10T 1 2 » 4< , D35X1  2«  4<»  D40U2,4C,  1)45X12»4<,  F'.X4<,  Fl0t4<,  Fl'>X4<,  f ft  10 
220X4C,  F25X4<«  F30X4C,  L5U2.4C,  L10X12.4<,  LlbT12.4<.  L20U2,4<,  ft  ll 
'L25X  12, 4<,  L JOt  l 2»  4<  , L J5tl2 ,4c  , 1.40X12, 4<,  L45tl2,4<  G 12 

RETURN  C.  13 

ENIRY  CMC  Ml  ft  14 

C Z.LRU  I ME  1NDICA1ED  ARRAYS.  ***  G L5 

DO  l K01.4  ft  16 

r5XK<0f  lOXKS#I  l5tK<0F2OXK<#|-25XR<*F3OXK<»O.O  ft  W 

DO  l JIM.'  ft  18 

06X J *Ks«l)l  OXJ  , K<0D  1 5t,l , KS #1)2 07,1 * Ke0U25  T.J  ,K4#l)30X  J ,K<  0D35tJ  »Kell)40f  ft  19 
l J,K<#n46tJ,K<»0.0  ft  »0 

L4X.I  »N\0L  lOiJ#  N<»01  ISXJ.KOl  2 At. I,  Kell  25XJ.KOI  U'*J,K<0L  TSXJ.Ktll  40X  ft  21 
l J,Ke#l  4S'*J,K<0O  ft  22 

4 T 
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J CONI INUE 

RETURN 

ENTRY  CATCH2 
IF  XLTSI.E0.5<  L5XL  >K<(  1 
05XL,K«f05XL,K<f.l 
IF  XLTST. LI. 10<  GO  TO  2 
- -IF  XLTSr.OO.lOC  .LIQXL.KOl. 
0l0XL,K<#010XL,K<£t 
IF  XLTST.LT. 15<  GO  TO  2 
IF  XLTST .CO. IS<  L15TL,K<#1 
UlSXL,K<it015XL,K<£l 
IF  tLTST.LT. 20<  GO  TO  2 

l£.XLLST.Efl.2Q.<  L2GXL.  R<«1.. 

D20XL,K<#020XL,K<U 
IF  XLTST. LT.2S<  GO  TO  2 
IF  XLTST .KQ.25<  L25IL,K<#1 
l)25XL»K<#025XL,K<f.  1 
IF  XLTST. LT.30<  GO  TO  2 

. IF  XLTST. EQ.3IX.  L30XL.K<il_ 

030X L»  K<#D30XL»  K<£ 1 
IF  XLTST. LT.35<  GO  TO  2 
IF  XLTST. F0.3S<  L35U»K<#1 
035XL,K<#U35XL,K<£1 
IF  XLTST. LT.40<  GO  TO  2 

- - if.  J.LIST..L0.40C  .L40XL,K<£1.. 

>40XL,K<tf(J40XL,K<f.l 
IF  XLTST. U.4S<  GO  TO  2 
IF  XLTST. EQ.45<  L45XL.KOII 
045  XL  * K< #045  XL  • K<£ 1 

2 CONTINUE 

. . .RETURN,  .... 

ENTRY  GILCAT3 

C GATHER  STATISTICS  ON  ONE-SLCONO 

C PROllABIL  IT  I6S.  *** 

K5#K10*0 

M54KK-4 

MIQKKKt.Q  . . _ 

00  3 Jrfl.12 
00  3 KK 1*2 
K5#K5C05XJ,K< 

K10*K10£010XJ,KC 

3 CONTINUE 

— S5*xxs*a.a</x24..*MSc  _ _ 

S10#XK10*1.0</X24.*M10< 

PRINT  6 

PRINT  7,  SStSlO 
LL5f LL10A0 
00  5 Kft 1 *2 

, UO — . 

4 l#I£l 

1 IF  X I .GT . 1 2<  GO  TO  5 

LL5XLL5SL5XI ,K< 

LLI0DLL1 OfiLIOTI «K< 

GO  TO  4 


G 23.. 

G 24 
G 25 
G 26 
G 27 
G 28 

. G — 29. 

G 30 
G 31 
G 32 
G 33 
G 34 

G 36 
G 37 
G 38 
G 39 
G 40 

G 42 
G 43 
G 44 
G 45 
G 46 

- . . — ..a.  .4;.-. 

G 48 
G 49 
G 50 
G 51 
G 52 

- G.  53 

G 54 

ANT  TWO- SECOND  LONG  CFLOS  G 55 

G 56 
G 57 
G 58 

- - - - ...  - . . _G~59.  - 

G 60 
G 61 
G 62 
G 63 
G 64 

G 66 
G 67 
G 68 
C.  69 
G 70 

... . C — 7-1— 

G 72 
0 73 

G 74 
G 75 
G 76 
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5 U1NUNUE  . .G.  7Z 

PR1NT  8 G 78 

_£K1M.  9..  LLZ.LLLQ G LS_ 

RETURN'  G 80 

.G  .8.1. 
G 82 

. G.  83. 

FORHAT  tlOX  • 54HRK0BAB1  LI  TY  OF  A CFLOS  FOR  AT  LEAST  THE  TIME  INDICA  G 84 

L!JiD< 

7 FORMAT  tl0X,5H.GE.t,F8.3,10X,5H.GE.2,F8.3<  G 86 

8 FORHAT  tlOX,74HMINlMUM  FREQUENCY  OF  OCCURRENCE.  OF  A CFLOS  FOR  AT  L G 87 

1 EAST  THE  TIME  INOICATEIX  G 88 

9 FORMAT  m.lOHONE  SECOND , 19, 3H/24.9X. IIHTWO  SECONDS, 19,3H/24<  G 89 

END  G 90- 
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LIST  OF  TECHNICAL  NOTES 

Number  Title 

74-1  Atmospheric  Moisture  Parameterization  (AD-784814) 

74-2  Development  of  a Gridded  Data  Base  ( ) 

(Publication  delayed) 

74-3  A Precipitating  Convective  Cloud  Model  (ADA-002117) 

74- 4  A Synoptic-Scale  Model  for  Simulating  Condensed  Atmospheric 

Moisture  (ADA-002118) 

75- 1  Estimated  Improvement  in  Forecasts  of  the  SANBAR  Hurricane 

Model  Using  the  Airborne  Weather  Reconnaissance  System 
(ADA -004097) 

75-2  Spring  Weather  Patterns  of  the  Western  United  States 
(Reprints)  (ADA-OO669I) 

75-3  Summer  Weather  Patterns  of  the  Western  United  States 
(Reprints)  (ADA-OO986O) 

75-4  Autumn  Weather  Patterns  of  the  Western  United  States 
(Reprints)  (ADA-013801) 

75- 6  Winter  Weather  Patterns  of  the  Western  United  States 

(Reprints)  (Publication  delayed) 

76- 1  Listing  of  Seminars  Available  at  Hq  AWS,  AWS  Wings,  and 

AFGWC  (Publication  delayed) 

76-2  Some  Aspects  of  Estimating  the  Probability  of  Cloud-Free 
Lines-of-Sight  in  Dynamic  Situations  ( ) 
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Date 
Jan  74 
Apr  74 

May  74 
Jun  74 

Jan  75 

Mar  75 
May  75 
Jul  75 
Sep  75 
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